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Proceedings of the American Society of Civil Engineers 


SELECTION OF DESIGN WAVE FOR OFFSHORE STRUCTURES®# 


Charles L. Bretschneider,* A.M. ASCE 
(Proc. Paper 1568) 


ABSTRACT 


Methods are presented for selecting for offshore structures the design 
wave, crest elevation, wave length and the maximum horizontal particle 
velocity at the crest.- The problem is somewhat different in regards to 
coastal structures—sea walls, rubble mound breakwaters, levees, etc., and 
will not be discussed in this paper. Essentially, this paper is the third of a 
series intended to supplement the information presented in an earlier report 
on Surface Waves and Offshore Structures by Reid and Bretschneider.(1) The 
first paper “Hurricane Design Wave Practices” by Bretschneider(2) presents 
methods for computing the significant and maximum waves for a hurricane 
moving from deep water over the Continental Shelf to the coast. The second 
paper “Approximate Response of Water Level on a Sloping Shelf to a Wind 
Fetch which Moves Towards Shore” by Reid(3) presents methods for obtain- 
ing hurricane storm surge. 

Wave data from various sources, supplemented by a limited amount of 
wave theory required for the development of generalized graphs are given in 
Appendix I. This information is used to revise some of the original graphs of 
reference (1). The wave data include some from the large tank at the Beach 
Erosion Board; some hurricane wave data from Lake Okeechobee, Florida;(4) 
model data by Danei(5) of France; that summarized by Munk;(6) and that re- 
ported by Mason(7) from the small tank at the Beach Erosion Board. 

Much of the material in the present paper was originally initiated by 
Professor Robert O. Reid of Texas A. & M. as presented in reference (1), and 
therefore, should be credited accordingly. 


Note: Discussion open until August 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1568 is 
part of the copyrighted Journal of the Waterways and Harbors Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. WW 2, March, 1958. 

a. Presented at ASCE Convention, New York, N. Y., October, 1957. 

* Hydraulic Engr. (Research), Beach Erosion Board, Corps of Engrs., U. S. 
Dept. of the Army, Washington, D. C. 
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INTRODUCTION 


Methods of predicting maximum waves for hurricane design wave prob- 
lems have been presented in reference (1). In regards to both the design of 
offshore structures and coastal structures, it is important to realize a design 
wave height, period, crest elevation and wave length. In many cases there 
may be a number of waves for selection. For the design of an offshore plat - 
form supported by piling, it is necessary that the underside of the deck clear 
the crest of the wave. The wave length and surface profile are important 
when computing total forces and overturning moments induced by the wave 
on a series of supporting piles. The maximum horizontal particle velocity at 
the crest gives a check against the wave theory used to compute the internal 
velocity field. 

The selection of the design wave consists of the following steps: 


I. Calibration 
(a) Storm Surge Formulae 
(b) Wave Forecasting Formulae and Techniques 


Il. Frequency Study 
(a) Wave Statistics 
(b) Hurricane Frequency Study 


Il. Predictions for Design 
(a) Storm Surge 
(b) Significant Wave Height and Period 
(c) Maximum Wave Height and Corresponding Wave Period 
(d) Crest Elevation of Maximum Wave 
(e) Wave Length of Maximum Wave 


The above steps are discussed below: 


Calibration 


It is necessary to calibrate storm surge prediction formulae and techniques 
to take into account deviations between observed and predicted high water 
levels. This is particularly evident in and near coastal areas where irregu- 
larities in both bottom topography and coastal features play important roles 
in modification of storm surges. 

Methods of predicting storm surges on the open coast are given by Reid(3) 
and Freeman et al, (8) among others. In using any method, storm surge hydro- 
graphs should be used where available for calibration, such that when applied 
to the design storm, a reasonable estimate of storm surge is obtained. This 
requires the analysis of all past hurricanes, where such data are available. 

In those cases where no data are available, a nearby and similar location 
where data are available may be used. As a first approximation to the solu- 
tion of the storm surge one may use directly the method proposed by Reid. (3) 

Calibration of wave forecasting formulae and techniques is also necessary 
to take into account variabilities in individual techniques and locations. A 
number of methods are available for predicting waves both in deep and shal- 
low water. In the case of deep water there are four methods generally used: 
Sverdrup-Munk-Bretschneider method, (9) Pierson-Neumann-James 
method,(10) Darbyshire method(11) and Suthon-Braceline method.(12) All the 


_ region from which the principal data were obtained. Perhaps, even better 
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above have incorporated therein inherent characteristics associated with both 
the procurement of data from various locations and also the method of analy- 
sis and interpretation of these data by the individual contributors. A very 
objective verification study of the above methods was made by Roll, (13) and 
the general conclusion was that each method works better for the particular 


verification might have been made provided the individual contributors made 
the forecast, each using their own methods and techniques. The graphical 
approach method, utilized by Wilson, (14) should be a useful tool for calibra- 
tion purposes. 

For wind waves in shallow water, there exist two methods: Thijsse(15) 
and Bretschneider, (16) each giving reasonably consistent results for shallow 
water of constant depth. The method of Bretschneider, (7) extended to the 
Continental Shelf, is calibrated, based on hurricane wind and wave data from 
Lake Okeechobee, Florida, (4) and some ordinary wind wave data from the 
Gulf of Mexico.(18) Although this method is used to predict hurricane waves 
over the Continental Shelf, improvements in the techniques can be expected 
when hurricane wave data become available from the Continental Shelf for a 
refined calibration. An example of computations illustrating this method is 
given in reference 2. 


Frequency Study 


Much data on wave statistics for various locations along the United States 
shores have been compiled by use of wave hindcasting techniques. The data, 
summarized in various Beach Erosion Board Technical Memoranda, are 
given in terms of the significant height and period. The importance of these 
wave statistics can be considered, where hurricanes are not a menace, for 
certain design problems other than offshore structures, and for planning 
marine operations. It must be emphasized that the maximum wave may be as 
much as 1.87 times the significant height. 

To date, little information is available on the frequency of hurricane waves, 
as far as actual observations are concerned. However, hurricane wave sta- 
tistics for deep water of the Gulf of Mexico have recently been compiled by 
Wilson(19) using the graphical approach to wave forecasting.(14) For most 
design problems these data must be brought across the Continental Shelf tak- 
ing bottom friction and refraction into account as further generation takes 
place. This report might be a source for obtaining the design wave for loca- 
tions in the Gulf of Mexico. 

The selection of the design wave can be obtained by use of an alternate 
method, wherein the design storm is first selected. This method consists of 
a frequency study of all hurricanes and other major storms passing a given 
length of coastline. Generally a 300 to 400-mile section of coastline is se- 
lected in order to obtain a sufficient number of hurricanes for a frequency 
study. The frequency study is then normalized in terms of 100 miles of coast- 
line. For example, consider the hurricanes of record which entered the 
United States coast between Jacksonville, Florida and Cape Hatteras, North 
Carolina. The distance between Jacksonville and Cape Hatteras is close to 
430 nautical miles. Between 1893 and 1955, 17 hurricanes entered this part 
of the coast. Insufficient data are available for hurricanes occurring before 
1893. These 17 hurricanes are summarized in Table A in descending order of 


x 
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TABLE A 


Frequency of 17 Hurricanes Between Cape Hatteras and Jacksonville 
(1893 = 1955) 


Date — yee of RAP (percent ) 
1 9/31/99 110 1.22 134.2 2.94 
2 10/15/54 36 2.31 83.2 8.82 
3 9/19/55 50 1.46 73.0 14.71 
h 8/12/55 LS 1.0 63.0 20.59 
5 12/ 2/25 1.02 55.1 
6 8/27/93 34 1.6 50.8 32.35 
7 9/13/13 39 1.16 45.2 38.23 
8 9/ 2/98 32 1.40 4.8 4.12 
9 8/17/55 LS 0.82 36.9 50.00 
10 9/17/06 37 0.99 36.6 55.88 
11 9/13/93 23 1.59 36.5 61.76 
12 8/11/40 27 1.19 32.1 67.64 
13 8/28/11 27 1.05 28.) 73.53 
14 8/14/53 30 0.79 23.7 79441 
15 10/15/47 12 1.38 16.6 85.29 
16 8/11/14 19.5 0.76 14.8 91.17 
17 8/30/52 20 0.95 ll.) 97.05 


(Values R and AP were furnished by Hydrometeorological Section, 
U. S. Weather Bureau 


} 
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the Energy Index, RAP. The frequency of occurrence of the above hurricanes 
is obtained by use of the Beard’s formula.(2) 


¢ = 100 (N 
M 


f = percent of occurrence (frequency) 


» where 


N = the order of RAP (1 to 17) 


M 


the total number of hurricanes, 17 


The information from Table A is plotted on the log-probability graph, 
Figure 1, and it is seen that the intensity along this part of the coast is slight- 
ly greater than for the United States as a whole (as given by Reid; see 
Figure 7, reference 2). 


Since 17 hurricanes entered a 430 mile section of coast in 62 years, one 
62 430 


obtains T7 * Too = 15.7 years as a hurricane recurrence interval, which 
means on the long run average that the center of a hurricane will enter a 100 
mile section of the east coast between Cape Hatteras and Jacksonville once in 


15.7 years. The recurrence interval for any particular strength of storm for 
Figure 1 is obtained from 


I= 15.7 100, or f = 1570 
f I 


Figure 1 can be used to obtain the Energy Index RAP for a particular re- 
currence interval, say once in 100 years, whence f = 1570/100 = 15.7 percent, 
and from Figure 1 RAP = 69. From Table I it is seen that the radius of maxi- 
mum wind will be between 45 and 50 nautical miles and AP between 1.46 and 
1.40 inches of mercury. Further investigations reveal that such a hurricane 
may have a forward speed VF between 7 and 14 knots (Table I, reference 2). 
Hence the design storm is selected for a recurrence interval of once in 100 
years, per one hundred miles of coast: 


R = 48 nautical miles 
AP = 1.4 inches of mercury 
RAP = 69 


Vr = 10.5 knots 


Design Wave Predictions 


Once the design storm has been determined and the prediction formulae 
calibrated, computations are made for storm surge and significant waves. 
Methods used are given in the references. An example is worked out in 
reference 2 based on Hurricane Hazel of October 15, 1954. For this example 
RAP = 83.2, corresponding to f = 9.0 percent from Figure 1. The recurrence 
interval in years per 100 miles of coast is I = 1570/9 = 175 years. In this 
example it was found for a mean low water depth, do = 74 feet, the total depth 
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with surge applied d = 85.3, and Hmax = 39.6 feet and T = 11.63 seconds. 
The crest elevation above d = 85.3 feet is obtained by use of Figure 2.* 
Compute H/2 = .293 and d/72 = .63, and read from Figure 2, 1 0/y = .69. 


Hence 1), = (39.6) (.69) = 27.3 feet, and the crest elevation will be 85.3 + 27.3 
= 112.6 feet above sea bottom. 


The wave length is obtained from 


L = gt? tanh , where 
“A La 


La is the Airy Wave Length and L/La is a correction factor for wave steep- 
ness and is obtained by use of Figure 3. The deep water Airy Wave Length is 


Lo = 5.12 T? = 693 feet. Compute d/y, = a3 = .123 and from Tables** by 
L L 
Wiegel, (21) tanh = .7656 (or from Figure 4 hn -77). From 
La Lo Lo 


Figure 3 for 9/2 = .63 and ¥/p2 = .293 read “/L, = 1.06. 


Thus L = 1.06 (693)(.77) = 565 feet. The wave celerity c = ©/ 
= 48.5 feet/second. 
The maximum menimane particle velocity at the crest can be obtained 
from Figure 5. Compute, (64.4427.3) 


= 565/ 


T 11.63 


= .75. From Figure 5 read 


uu = 0.5. Thus u = 0.5 (48) = 24 feet/second. 
max/C max 
It is interesting to compare the above values with those obtained by use of 
the Airy(22, 23) wave Theory. No= a = (39.6) jo = 19.8 feet, or 7.5 feet lower 


than that given for this example. The Airy Wave Length is La = 693 (.7656) = 
531 feet or 32 feet shorter. The wave speed is C = L/>p = 45.6 feet/seconds. 


The maximum horizontal particle velocity is given by Rea * ou = 13.8 feet/ 
second. The Airy theory is usually extended to include the finite height, from 
which — is given by 


= 17.0 feet/second 


which is 40 percent less than that obtained by use of Bernoulli’s equation 
(Figure 5). Since the force on a pile is proportional to the velocity squared 
one obtains (24 44)? = 2.0 as the ratio of maximum crest force obtained by 

the two methods, assuming the same drag coefficient in either case. 

Thus if the Airy Wave Theory is used for design conditions, an added safety 


* See Appendix I for development of Figures 2, 3 and 5. 
The solution of La/t, tanh i is given in Figure 4 in terms of d /p2. 
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Umax 


FIGURE 5. MAXIMUM HORIZONTAL PARTICLE 
VELOCITY AT CREST 


factor is mandatory for both the crest elevations and the velocity and acceler- F 


ation fields. Use of the Airy Theory in the above sense is not a recommended 
practice. 


Application to Hurricane Wind Waves 


The preceding discussion is applicable to regular and simple waves, unaf- 
fected by strong winds. Figure 2 (and also 3) does not take into account any - 
change in wave characteristics due to strong winds. This can best be studied 
by use of emperical data. Hurricane wave data from Lake Okeechobee, ; 
Florida,(4) have been analyzed for H, T, d, and No (Table I-3). Figure 2, with 
the aid of H/y2 and d/p2 computed from Table I-3 was used to obtain the pre- 
dicted value of 7), also presented in Table I-3. Figure 6 is a comparison be- 
tween observed 1), and predicted 7,4. Much of the scatter is due to the ina- 
bility of obtaining exactly, H, T, and d, and long run averages are more 
significant. Even under controlled laboratory conditions one experiences 
considerable scatter of individual data. The regression line of observed 15 
to predicted 7, has an intercept at .052 and a slope of 1.045 as determined by 4 
the least squares techniques. If it is assumed that Figure 2 is correct for a 
waves unaffected by wind, then the effect of hurricane winds is to increase on \ 
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Regression Line 


DATA FROM LAKE OKEECHOBEE 


Loke Station 12 
Lake Station 14 
Lake Station 


Observed crest elevation, feet, [7] 


Equation of Regression Line 
052 + 1.045 [%] 
ob pr 


Correlation coefficient, r = .958 


° 1 2 3 4 5 6 
Predicted crest elevation, feet [No] 
pr 
FIGURE 6. COMPARISON BETWEEN OBSERVED AND PREDICTED 
CREST ELEVATION ABOVE STILL WATER 


(Predictions by use of Figure 2 ) 


long run averages the crest elevation above still water by about 5 percent. 
Because of the scatter of data, perhaps a 10 percent increase in some cases 
is justified for the design wave. This, of course, needs further study, using 
possible wave data of the Gulf of Mexico. 


SUMMARY AND CONCLUSIONS 


This paper represents the latest revision for obtaining for the design wave 
height, crest elevation, wave length, and maximum horizontal particle velocity 
at the crest of the design wave, and is intended to supplement the information 
presented in reference 1, as well as to be a continuation of reference 2. 

Two generalized graphs are prepared using wave theory, Bernoulli’s equa- 
tion and wave data. Figure 2 is for obtaining the crest elevation and Figure 3 
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for obtaining the wave length of the design wave. A fourth graph, Figure 5, 
can be used in conjunction with Figures 2 and 3 for predicting the maximum 
horizontal particle velocity at the crest. The effect of hurricane winds has 
been considered. Based on Lake Okeechobee hurricane wave data it was 
found that the crest elevations were about 5 percent greater than that predict- 
ed from Figure 2. The scatter of data shows that perhaps a 10 percent in- 
crease may be justified, but this needs further investigation of larger hurri- 
cane waves, such as might be obtained in the Gulf of Mexico. 
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LIST OF SYMBOLS 


A wave amplitude above still water, feet 


Ao H/o, half the wave height, feet 


b used as a subscript to denote a breaking wave 

Cc wave celerity feet/second 

Ca wave celerity of the Airy wave, feet/second Bi 
Ch wave celerity of breaking waves, feet/second 

d depth of water, feet 

f frequency of occurrence, percent 

g acceleration of gravity, 32.16 feet /second2 


trough elevation above bottom, feet 
wave height, feet 


Hp breaking wave height, feet 


I recurrence interval, years 
wave number, k = 
L wave length, feet 


Lo deep water Airy wave length, feet 


La Airy wave length, feet 


R radius of maximum wind, nautical miles 
yy wave period, second 7 
Us horizontal particle velocity along the free surface of gravity wave 


(stationary coordinate system), feet/second 


4 h 
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Up horizontal particle velocity at crest (@ = 0), feet/second 

Umax maximum horizontal particle velocity at the crest, feet/second 

VA absolute particle velocity, feet/second 

Ve forward speed of hurricane, knots 

Ws vertical particle velocity, feet/second 


Wm particle velocity at still water elevation at the free surface, feet/ 
second 


3.1416 
angular frequency, radians/second 
elevation of free surface of wave measured from still water, feet 
No crest elevation above still water, feet 
Nb crest elevation of breaking wave above still water, feet 


phase position, degrees 
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APPENDIX I 


The following material is limited to the development and application of 
Figures 2, 3, and 5. This requires only a brief summary of wave theory, the 
use of Bernoulli’s equation, and the source and nature of the wave data. The 
section on wave theory does not consider the internal velocity and accelera- 
tion fields usually reserved for papers on wave forces. If necessary, such 
information may be obtained from the cited references and numerous other 
sources. 


Brief Summary of Wave Theory 
(As applicable to Figures 2, 3, and 5) 


Airy Wave Theory 


The Airy Wave Theor (7, 8) is applicable to very low steepness waves in 
deep or shallow water (#/, and H/q small). The surface profile is given by 


1-3 cos @ , where 


(1) 
@ = kx -c t is the phase position 
= wave height, vertical distance between crest 


and trough 


wave number 
angular frequency 
horizontal distance from crest 


time position from crest 


wave period 


Airy wave length 


The Airy wave length is given by 


La er? tanh 2rd 
La 


er 
T 


DESIGN WAVE 
g = acceleration of gravity 
d = mean water depth 


The wave celerity, Ca is given by 


La 3 
(3) 


The Airy theory is used for the lower boundary limit (Figures 2 and 3) in 
2. the field of wave theories. 


Solitary Wave Theory (for breaking waves) 


The solitary theory(6) for non-breaking waves will not be discussed. The 
highest wave in shallow water is governed by the solitary wave theory(6) and 
is a wave of elevation entirely. The height of the breaking wave is obtained 
from Hp = 0.78d. 

The celerity C), is given by 


Ch = VY 2g Hy (4) 


The crest of the breaking wave is cusp shape with angle of 120° and the 
crest elevation No/y = 1.0. 


This theory is used for the upper boundary limit for dy > 0 


Modified Solitary Theory 


The modified solitary theory discussed in great detail by Munk, (6) is ap- 
plicable to breaking waves in very shallow water. The breaking wave height 


is obtained from Hp = 0.78h, where, h = elevation of trough above bottom, 
given by 


(5) 


. The wave celerity of the breaking wave is given approximately by equation 
4, and more exactly by equation 18. The wave celerity for a non-breaking 
solitary wave is given by C = /g(H + h), and applies only for moderate or low 
values of 8/H. The wave length is obtained from L = CT. 

The elevation of crest above still water is obtained from 


(a - h) (6) 


( H ) 


This theory is used as the upper limit in very shallow water, small dy, . 


Michell - Havelock Theory 


The Michell-Havelock theory(24) is applicable to breaking waves in deep 
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water. The breaking wave height is given by, Hp = L/7, where L is the wave 
length, given by 
L= 1.20 (7) 
27 


The crest elevation above still water is obtained from NO/y = 0.676. 
This theory is used for the upper limit in deep water, d/y = 0.5. 
Stokes Theory (second approximation)(25) 


This is a second order theory, one order above the Airy theory, and ap- 
plies for waves of small finite steepness. The wave length is given by 


2 
L A (tanh) kd) 1+ (yf (F) » where 


d 5 + 2 cosh 2kd + 2 [cosh [2ka | |? 
f, (F) 
L 
8 [sinh kd | 
The surface profile is given by 


Wy =i cos fo (9/;) cos 20 , where (9) 


2 


(cosh kd) (2 + cosh 2kd) 
[sinh ka] 


The elevation of crest 1), above still water is obtained from 9 for @ = 0. 
This theory is used next to the boundary of the Airy theory. 


Stokes Theory (third approximation for deep water) 


The third approximation for deep water has a surface profile given by 


7 wh 2 
Wy = cos @ + cos 20 + cos 3@ where (10) 


A is related to Ag = 4 according to Keulegan(26) 


wA 2 11 


and the celerity is given by 


(12) 


This theory is used as the boundary for deep water, for waves of finite 
height, but is adjusted slightly as H/; approached I/n, 


| 
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The solution of equations 11 and 12 is given in Table I-1. Equation 10 for 
N o/y breaks down sooner than 12, and is adjusted by use of wave data. 
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TABLE I-1 
Solution of Stokes Third Approximation i, 
for Wave Length in Deep Water ir 
1.027 0275 
20748 1.054 2397 
20925 1.083 
21079 1.112 
01217 1.141 
1.172 
The above theories are those used, together with wave data to prepare a - 
set of useful graphs for determining wave characteristics for most design 
waves. The Cnoidal wave theory might also be mentioned, since it represents a 
waves in the intermediate zone between Stokian and solitary waves. However, 
this theory is a little difficult to use in its present form, and instead is re- 


placed by wave data as far as Figures 2 and 3 are concerned. 


Application of Bernoulli’s Equation 


Wave data with the aid of Bernoulli’s equation can be used to obtain best fit 
numerical relationships for wave characteristics in the intermediate zone 
where no simple wave theory exists. The problem of wave motion can be re- 
duced to one of steady state by use of the Rayleigh principle wherein a steady 
current equal to the wave celerity is superimposed thereon in an opposite 
direction to wave motion. The free surface profile then becomes a stream- 
line. The pressure along this streamline is constant atmospheric, or in 
reference to the atmosphere is zero. Using the mean water depth as a 
reference and applying Bernoulli’s equation along the free surface of the now 
stationary wave, one obtains 


"| « « £g% , where (13) 


> 
| 
: 
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(u,-C) is the horizontal particle velocity along the free wave surface in the 
moving coordinate system, ft./sec. 
ug =the horizontal particle velocity along the free wave surface in the 
stationary coordinate system, ft./sec. 
C =the wave speed, ft./sec. 
Wg = the vertical particle velocity at the free surface, ft./sec. 
n = elevation of free surface measured from still water (positive above 
and negative below still water), feet. 
g = acceleration of gravity 32.2 ft./sec.2 
K_ =the constant of Bernoulli’s equation 
The absolute velocity Va, is obtained from 


= (ug - Cc)? + We? (14) 


The constant K can be evaluated at any number of known positions of the 
free surface. For the purpose of this discussion, the constant will be evalu- 
ated at still water, where @ = 0. At this point ug = 0, but wg has a finite value 
which will be termed wy, not necessarily the maximum vertical particle 
velocity. Hence the constant becomes 


2 
K=el¢ (15) 


(16) 


The above expression relates the surface profile to the wave celerity, hori- 
zontal particle velocity ug and wm. 

The crest elevation 1, above still water can be determined, since at 8 = 0, 
Ug = Up and wg = 0, whence 


22 No 2 + (17) 


m 


Breaking Waves 


For the special case of breaking waves (deep or shallow) 


Uy = C=C 


2 
= 26 Ym 


Equation 18 can be written in the following form 


Thus 

2 2 2 

thus 
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Cp 
H/ 2 
T 


Where La is the Airy wave length given by equation 2. 
For the special case of the solitary breaker in shallow water, 


Mo 1.0, Wine oO, tanh — one 


and H/q = .78, whence from equation 19, L/L, = 1,25. 


For the special case of the highest wave in deep water, according to the 
Michell Theory, 1b/y = .676, //L,a = 1.20 and H/y = 1/7 and since tanh 


w 
ana > 1.0, whence from equation 19,(<—)* = .0113. 
b 


Ww 
If it is assumed that ra 2 varies from 0 is shallow water to a maximum 
b 


of .0113 in deep water, then its contribution in equation 19 is negligible for all 
breaking waves. 
Thus 


Cp = Jem, for all breaking waves. (20) 


Maximum Horizontal Particle Velocity 


The maximum horizontal particle velocity, u,,,,, occurs at the crest posi- 
tion @ = 0 (Umax = Yo, and N =o). Umax Can be investigated by use of 
Bernoulli’s equation, provided one has previously determined 7, and L, say 
from Figures 2 and 3 of the text, respectively. Solving equation 17 for 
Umax/c one obtains 


21 
Umax = 1 + (21) 


ce 


It was shown that for all breaking waves deep or shallow water (wm/o)” 


can be neglected. This is also probably true for all non-breaking waves, 
whence, 


Mmax 21 (22) 
V 


Umax Can never be greater than given by the above approximation. 
Equation 22 can be investigated for the Airy Wave Theory, H/L very small. 


If a is very small, 1, = , and gH/o2 << 1.0. Hence an expansion of equa- 


tion 22 by the Binominal theorm leads to 


ag (23) 


q 
- 
q 
; 
oH 
cu 
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which is the same as that predicted from the Airy theory for H/; very small. 
consider also equation 21, for the Airy Wave Theory, (Wm)2 = (Umax)* and 

one obtains exactly equation 23. Thus equation 22 is a very good approxima- 
tion for obtaining the maximum horizontal particle velocity at the crest. The 
solution of equation 22 is given in Figure 5 of the text. 


Source and Nature of Wave Data 


Wave data used in this paper consists of both breaking and non-breaking 
waves. Breaking waves from various sources have _ summarized by 
Munk(6) in terms of Hb/,? and dye where Lo = ‘ T*. This information 
in terms of Hp/p2 and d/p2 is given below in Table I-2. 


TABLE I[-2 
(After Munk) 


Summary of Observations of Depth of Breaking Location 


~ Number of 


Y 2 d 2 Symbol 
Observations /t 2 In 
Which Were feet/sec* feet/sec Fig. 1 


Averaged 


Estero Bay 2hh 2026 O34 


O 
SIO Daily 238 20) 0053 A 
SIO Leica 56 2007 A 
WwHOI 17 .072 .079 
Berkeley Tank 2 0123 2157 
B.E.B. Tank 27 +205 2207 
Lake Superior 134 2358 0558 


The above information is presented in Figure I-1, Breaking Index Curve, 
for H/yp2 versus A/p2, At the far right of Figure I-1 are two points; (symbol 
¢), wind obtained in the Berkeley tank by Bretechasidet and Rice.(28) 

The other symbols (*) are data from Danel,(5) originally given in terms of 
H/, and 4/;, but transformed approximately to H/p2 and 4/2. The lower 
portion of the curve is governed by the Solitary theory (H = .78d) and the deep 
water portion by the Michell theory (H/; = 1/4, or H/72 = .88). The transi- 
tion is governed by wave observations, imposing Bernoulli’s equation and the 
other figures to be discussed. 

Figure I-2 is the Breaking Index Curve in terms of 4/, versus 4/L and 
is similar to that of Figure I-1. The dashed curve is the relationship given 5 
by Danel(5) and corresponds to H/L = 1/q tanh ona . At about d/y = 0.1 the 
curve by Danel intersects and crosses the solitary limit H = .78d. The solid 
curve is a continuous curve from 4/; = 0 to 4/; = deep water, and is based 
on the imposing solution of Bernoulli’s equation using Figure I-1 and others 
following. 


Figure I-3 is a plot of the breaking crest elevation above bottom Yo/p* 


versus 4/52, and the data is that obtained from Lake Okeechobee, Florida. 
This curve is obtained by use of Figure I-1 and Figure 2 of the text, where 
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Yo /p2 = q/,.2 + No (H/-p2). The data consists of maximum waves of the 
H 


October 17-18, 1950 hurricane.(4) These waves are not necessarily breaking 
and there is no way of knowing which waves were breakers. Hence, the enve- 
lope curve and not a mean curve is used. The data are summarized in Table 
I-3, and not all the points are plotted in Figure I-3. 

In order to tie in the breaker index relationships, with the figures follow- 
ing, equation 19 was used. Figure I-4 is an imposing solution of this equation 
assuming Figures I-1 to be correct, and assuming (Wm/c,,)2* varies as 


shown from 0 at Hb 2 = 4/p2 = 0, to .0113 at deep water, H/p2 = .88. As- 
suming Figure I-2 to be correct permits determination of L/Ly as a function 


of either 4/T2 or H/72, and hence the determination of Mb/y for all breaking 


waves in the transition zone. Perhaps other representations of these figures 
could have been made, thus altering each of the corresponding figures. How- 
ever, there is one best solution which gives the most accurate overall rela- 
tionships. A trial and error process was used to obtain this best solution. In 
addition use was made of the non-breaking waves. 

Figures I-5, -6, and -7, plots of 7, versus H, are wave data for non- 
breaking waves obtained in the large tank of the Beach Erosion Board. The 
depth of the water was 15 feet. Three wave periods, T = 5.6, 7.87, and 11.33 
seconds, were tested for a wide range in wave heights. Each wave record 
was 15 minutes in duration and the highest and the lowest individual wave of 
each record was determined. These data are summarized in Table I-4 ac- 
cording to the highest and the lowest wave of each record. The termination 
point of the curves of each of the above figures is obtained by the Breaking 
Index Curve of Figure I-1. These data are truly in the intermediate zone of 
and H/p2: 

Figure I-8 is a plot of data of Toya versus H/y2 for deep water waves 


obtained in the small tank of the Beach Erosion Board, and reported by 
Mason.(7) Since these are deep water data the variable d/,,2 is eliminated. 


However, the wave period T was different under each test run, hence the plot 
is given in terms of No/H versus H/T2, The breaking point corresponds to 
the Michell Theory, 10/y = .676 and H/p2 = .88. These data are summarized 
in Table I-5. 

Using cross plots of Figures I-5, -6, -7, and -8 and the previous figures 
discussed, it is possible to construct the generalized graph of Ng, _, Figure 2 
of the text. The graph of NO/y versus d/p2 with isolines of H/p2, are tied 


together as best fit relationships using wave theory, Bernoulli’s equation and 


wave data. Table I-6 represents the final selected values of the breaking 
wave characteristics. 


Figure 3 of the text, L/ta, needs some discussion. The envelope curve of 
i is related to the Breaker Index Curve, Figure I-1. Bernoulli’s equation 
Figure I-4, and the envelope curve of No/y of Figure 2 of the text. These are 


* It would have made very little difference in the final results if it had been 
2 = > 
assumed (Wm/c,) 0 for all breaking waves. 
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TABLE I-S 
Summary 
Deep Water Waves (Beach Erosion Board Small Tank) 
Run H( feet) T(sec) H/p2 
3 e264 2848 367 2610 
1 359 «602 
Average 0379 0576 
18 2236 228 2610 
023k 2860 2316 2556 
2860 0315 0579 
13 e225 81.8 e213 587 
2 305 2988 e311 
23 e360 30 4 
0296 0988 302 
e257 0293 e540 
Average e323 2567 
19 2203 286 2589 
] 19 ol 2842 e541 
19 2192 2270 
178 e251 2539 
4 
Average 0275 e550 
2 1.003 e2h2 538 
lh ecat 969 487 
22 e170 e530 
1) 2218 969 0233 550 
1h ecll 969 2228 2509 
2 2220 1.003 563 
22 215! 2555 
5 ec07 2985 e2lh 
17 0192 2966 20k, e551 
Average e231] 2537 
0195 2985 2202 253 
17 e186 2966 2198 e506 
20 0176 0996 2178 
20 e161 2996 e162 526 
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terminating points. The starting point for low but finite steepness is given by 
Stokes’ theory equation 8. For waves of finite height in deep water, Stokes’ 
third approximation is used, equations 11 and 12. The values given in Table I 
are very slightly adjusted such that L/La = 1.20 at H/p2 = .88. These are 


used to obtain the isolines of H/p2 at the far right of Figure 3. To obtain the 
continuation of these isolines of H/,,2 to the breaking limit established along 
the envelope curve one could make use of Figure 2 and solving equation 16 for 
non-breaking waves. This would require the knowledge of the horizontal par- 
ticle velocity at the crest, since ug # C for non-breaking waves. At present a 
there is no suitable wave theory for obtaining up for large non-breaking 

waves. However, one could start with an assumed velocity field based on an 
approximation of wave theory, using Figures 2 and 3 as first approximations a 
and employ relaxation techniques to obtain the complete picture. Such a pro- 
cedure has been set up (Chappelear)(27) in case of the Solitary theory. The 

isolines of H/72 in Figure 3 have been assumed as shown, but have not been 

checked and adjusted according to relaxation techniques. It is believed that 

the assumed isolines of H/p2 of Figure 3 are sufficiently accurate for design 
purposes, but if greater accuracy is required one may refer to the work by 
Chappelear.(27) Hence the design wave length is given by 


Figure 5 is used to predict the maximum horizontal particle velocity at the 
crest, the accuracy of which depends on the accuracy of Figures 2 and 3. 
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SHIPBOARD HYDRAULIC BREAKWATER® 


R. A, Dilley! 
(Proc. Paper 1569) 


SYNOPSIS 


A model hydraulic breakwater was constructed on a 1:86.5 scale model 
Liberty Ship. The ship was moored in a model basin in seven inches of water 
and subjected to a series of different wave conditions. The mooring forces 
and wave attenuation were measured. 

Graphs showing the relationship between wave attenuation and jet-flow rate 
and the dependence of mooring force on wave height are included. It was 
found that the jet-flow rate to cause a given wave attenuation increased as the 
ratio of wave-length to water-depth increased. There was a critical wave 
height for each wave period at which the mooring line forces began to increase 
very rapidly with increasing wave height. 
The investigation was conducted in a 150 ft. x 64 ft. x 2-1/2 ft. model basin 


at the Engineering Field Station, University of California, Richmond, Cali- 
fornia. 


INTRODUCTION 


Many investigations have been made into the problem of providing a mobile 
breakwater which would be effective for calming normal ocean waves. The 

problem is complicated by the fact that in general the area to be calmed is in 
an area where the water is shallow and the wave action extends over the whole 
depth. Model studies and some prototype installations have shown that most 

of the proposed types of mobile breakwaters will reduce wave heights appreci- 
ably when the ratio of wave-length to water-depth is relatively small. 


Note: Discussion open until August 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1569 is 
part of the copyrighted Journal of the Waterways and Harbors Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84 No. WW 2, March, 1958. 

a. Submitted in partial fulfillment for the Master of Science in Engineering, 

University of California, Berkeley, August 1957. 


1. Junior Engr., Univ. of California, Inst. of Eng. Research, Berkeley, Calif, 
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However, the longer waves cannot be reduced in height without causing pro- 
hibitive forces on the mooring lines. gamer study has been conducted 
recently on hydraulic breakwaters. 2,3) Most of the work has been done with 
two-dimensional models. Very little has been done in three dimensions and 
hardly any work has been done on the forces on the structure due to the jet 
action., The few studies on three-dimensional models show that the waves 
tend to refract around the ends of the breakwater and augment behind it. 

One method of preventing this refraction would be to create currents in 
both directions. Instead of just directing jets into the train of waves, jets also 
would be directed in the opposite direction. This arrangement would have an- 
other advantage since it would eliminate the large thrust on the breakwater 
due to the jet action in only one direction. 


General Considerations 


A portable hydraulic breakwater would be very useful for protecting a 
section of ocean during an amphibious landing or some other temporary oper- 
ation. One way to have the hydraulic breakwater mobile would be to install it 
ona ship. Then it could be moved to any other location in the minimum of 
time. While in operation the ship would have to be moored so that it would 
remain in a specified position. 

In order to study the feasibility of a mobile hydraulic breakwater a model 
was constructed and then tested in a model basin (Fig. 1). A model of a Liber 
ty Ship hull was available for this purpose. The hydraulic breakwater 
equipment was installed on the model ship and a four mooring line system was 
designed to hold the ship in position while the jets were in operation. 

The problem of mooring a mobile hydraulic breakwater is more complicat- 
ed than that of mooring a standard ship due to the jet action. The compli- 
cations due to the jet action will be somewhat alleviated by having jets on both 
sides of the ship; however, the jets could affect the resonant period of oscil- 
lation of the ship model. Also it is likely that the jets will cause net forces 
on the ship due to slight variations in flow rate from one side to the other. 

The problem of mooring a floating structure such as a ship in the open sea 
where it is subjected to wave action is analogous to the problem of a mass 
restrained by springs and acted upon by some force which varies with time. 
However, due to the many degrees of freedom and the non-linearities of the 
mooring system no general theoretical solution has been formulated for the 
problem of a ship moored at any heading to the seas. Neither has there been 
any theoretical solutions to the relatively simpler problem of a ship in head 
or beam seas. The pioneer effort in a theoretical solution of this problem was 
made by Wilson, 1950.(6) Wilson developed a differential equation for the 
longitudinal motion of a ship moored alongside a dock and subjected to rela- 
tively long period harbor seiches. Wilson’s analysis is considerably simpler 
than that which would be necessary for a vessel moored in normal sea waves, 
since the particle motion of a seiche is more easily described than that of a 
wind wave or swell and also since he only considered longitudinal motion. A 
ship moored in a particular heading to the waves is not limited to motion in 
only one direction. However, certain relationships were developed by Wilson 
which yield some valuable qualitative results. As would be expected the reso- 
nant period of oscillation varies inversely as the initial tension in the mooring 
lines. The amplitude of the seiche and consequently of the vessel motion 
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affects the resonant period of oscillation when the tensions are high; the 
greater the seiche amplitude, the lower the resonant period of the system. 
Wilson and Abramson, 1955,(7) have used the Ritz method to treat the Wilson 
equation and have shown some results for ship displacement as a function of 
a non-linear natural period of oscillation function for various seiche ampli- 
tudes. Although the above-mentioned problem solution is a gross simplifi- 
cation of the problem at hand it does show the effects of mooring line tension 
and possibly wave amplitude. 

Several studies have been performed(1,4) with different models of existing 
ships with several different mooring configurations. These studies indicate 
that the qualitative results mentioned above hold for ships moored in normal 
sea waves. The natural periods of oscillation vary inversely as the tensions 
in the mooring lines and the resonant period of oscillation becomes lower at 
higher wave amplitudes. 

There are many variables to be considered in the problem of mooring the 
mobile hydraulic breakwater. For the purpose of this first investigation only 
the following variables were considered and all others were held constant: 


(1) Seaward wave height 

(2) Leeward wave height 

(3) Wave period 

(4) Jet-flow rate 

(5) Mooring line tension (during one special part of the tests) 


Laboratory Equipment 


The experiments were conducted in a 150 ft. x 64 ft. x 2-1/2 ft. model 
basin at the University of California Engineering Field Station, Richmond, 
California. This basin is equipped with a piston-type wave generator at one 
end and a wave absorbing beach at the other end. A structure supporting a 
photographic crane spans the width of the basin affording a working platform 
which could be positioned over any portion of the basin. Suspended from this 
platform and immediately surrounding the area in which the tests were per- 
formed was a plywood and canvas wind protector. An instrument shed at the 
side of the basin housed the Brush Universal Analyzers and Brush Oscillo- 
graph used as recording equipment and the wave generator controls (Fig. 1). 

The water surface time histories from which the wave height and wave 
period could be determined were recorded with parallel-wire-resistance 
gages. 5) There were connected in parallel with one leg of a Wheatstone 
Bridge which was connected to a Brush Universal Analyzer. 

The model ship was a 1:86.5 scale model of a Liberty Ship which was out- 
fitted with four mooring lines, four force meters and the pump, manifold and 
jets which comprised the hydraulic breakwater system. With this equipment 
installed on the ship it was made dynamically similar to a Liberty Ship in the 
loaded condition according to Froude model laws. The apparatus for de- 
termining the ship characteristics is shown in Fig. 2. The longitudinal di- 
mension of the center of gravity was determined by suspending the ship in an 
upright position at two points, one near each end, and determining one re- 
action force. Then knowing the displacement (the ship was weighed) the hori- 
zontal position of the center of gravity could be calculated by taking moments 
about one point of suspension, The vertical position of the center of gravity 
was determined in the same manner, only the ship was supported on knife 
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edges on its side and the reaction force at one point was measured. The draft 
is the distance from the baseline of the ship to the waterline at the center of 
gravity and was measured directly. The draft varied along the ship because 
the model in the water had a slight slope downward toward the stern. The 
natural period of roll was measured directly by rotating the model in the 
water and allowing it to oscillate and timing a number of oscillations. The 
roll period presented is for the ship without mooring restraints. The radius 
of gyration for pitch was determined by suspending the model as a bifilar 
pendulum and allowirz it to oscillate in a yawing motion. From the period of 
the oscillations one can determine the radius of gyration by the following 
equation? from (4): 


where 


Tp = period of oscillation of bifilar pendulum 
g = acceleration of gravity 

= radius of gyration 
a&vw are shown in Fig. 2a 


This is actually the radius of gyration in the yawing direction but the two radii 
of gyration are assumed to be the same. The other characteristics tabulated 
were determined by direct measurement. The data for the various charac- 
teristics of the ship model are shown in Table I. 

he mooring force meters were designed to serve two purposes (1) record 
the force in the mooring lines and (2) provide the proper mooring chain elas- 
ticity. The force meters were U shaped double cantilever type with one 
Baldwin-Lima-Hamilton SR-4, A-5, 120 ohm resistance type strain gage glued 
on each side of the two cantilever arms (Fig. 3b). The force was recorded by 
arranging the four resistance type strain gages in a Wheatstone Bridge so 
that the maximum cumulative strain could be recorded for a given force on 
the meter. The bridge circuit could be hooked directly to a Brush Universal 
Analyzer which served as a power supply and amplifier. The force meters 
were calibrated for deflection characteristics and force measuring sensitivity. 
In order to give the mooring line the correct stiffness the meter had to have 
a given spring constant. The spring constant could be adjusted by changing the 
length of the cantilever arms. This constant was checked by applying a series 
of known weights and recording the deflection of the force meter at the point 
of the mooring line attachment. When the spring constant had been adjusted 
to the correct value, the length of the arms was fixed securely. 

In order to determine the force on the force meter from oscillograph re- 
cordings it was necessary to have a correlation between force on the meter 
and the deflection of the recorder pen. This was determined by statically 
loading the meter with known weights and observing the pen deflection. 

The mooring lines were made of No. 18 steel jack chain which weighed 
0.019 pounds per foot. This corresponded to a prototype 3-3/4 inch dielock 
chain. The stiffness of the model chain exceeded the value required as de- 
rived from model laws so that the force meter could be adjusted to give the 
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TABLE I 


MODEL SHIP CHARAC TERISTICS 


Center of Gravity 
Vertical, distance above base line, ft. 
Longitudinal, distance aft of bow, fte 
Draft and Displacement 
Draft at center of gravity, fte 
Displacement, lb. 
Natural period of roll, sece 
tadius of gyration, transverse axis, pitch, 
Mooring line positions 
Bow to bow mooring hawse pipes, fte 
Bow to stern mooring hawse pipes, ft. 


Baseline to mooring hawse pipes, fte 
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necessary extra elasticity. This requirement governed the length of the canti- 
levers used in the force meters. The mooring lines were attached to the 
meters with a short section of 12 mil wire. These were inserted through 
hawse pipes made of one-eighth inch copper tubing inserted through holes 
drilled in the sides of the model ship a short distance above the water line and 
clamped to the force meter bar. The hawse pipes were drilled at 45° to the 
longitudinal axis of the ship and with a dip angle of 45° below horizontal (Fig. 
3). The angle which the chain made with the ship varied as the ship moved, 
but the meter was constructed so that the force on it always acted perpendicu- 
lar to the meter arm. The mooring lines were 42 inches long and had small 
flat plates at the anchor ends. These plates had slots in them to allow for ad- 
justing the initial tensions. The plates were fastened securely to the basin 
floor with Phillips shields and screws. 

The hydraulic breakwater system consisted of a small pump and direct 
current motor, a one-inch plastic tube manifold running the length of the ship 
and ten jets, one-quarter inch in diameter on each side of the model. The jets 
were connected to the manifold with flexible tygon tubing and were fastened to 
the ship with slotted plastic plates and screws so that they could be easily ad- 
justed in the vertical direction. The jets were set to operate at the waterline 
during all the tests. The discharge of the jets could be varied by changing 
the speed of the pump which was varied by changing the input voltage to the 
pump motor. The pump was calibrated so that the flow rate was known as a 
function of motor speed and of motor input voltage. The voltage was recorded 
as a measure of the flow rate during the runs since the voltage was easier to 
measure during a run than the motor speed. 

The pump used was a surplus aircraft fuel transfer pump. Flow rates from 
20 to 70 gallons per minute were attainable at heads from 5 to 25 feet of water. 
The pump was powered by a one-half horsepower direct current motor. The 
direct current power was supplied by a motor-generator set with field current 
control to supply the variable direct current voltage source. 


Experimental Procedure 


Quantitative measurements were made of the mooring line forces at the 
ship, the incident and transmitted wave profiles and the jet-flow rate. The 
surface time history recorder between the wave generator and the model was 
placed 21 feet from the model so that the disturbance due to the model had a 
negligible effect on the wave profile at the recorder. This insured that the 
surface time history measured was the true incident surface time history. 

The wave recorder located on the lee side of the breakwater was situated in 
what was determined by several preliminary runs to be the area of maximum 
wave attenuation 8 feet behind the model. 

Before the tests were begun several preliminary measurements and cali- 
brations were made. The ship model with the hydraulic breakwater on board 
was assembled and the ship characteristics measured. The pump and mooring 
force meters were calibrated. These procedures are described in the section 
entitled “Laboratory Equipment” above. 

Besides the quantities which were measured for each run there were sever- 
al quantities which were measured periodically throughout the runs. The 
natural periods of oscillation in the three modes of surge, sway and yaw were 
measured. This was done by displacing the ship model in a direction so as to 
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cause the desired motion. The motion induced forces in the mooring lines 
and made a force time history on the oscillograph from which the period of 
oscillation could be measured, 

When the model was set up for each series of tests the initial tension in 
each mooring line was adjusted to a predetermined value. At intervals during 
these tests these tensions were checked to insure that they were maintained 
correct. During the tests the procedure was as follows: 


1. The wave generator and recording oscillograph were started. 
The jets were started and the motor voltage recorded. (When taking 
runs with jet action.) 


3. The mode of oscillation experienced by the ship was observed. 


The test runs covered a wave period range of 0.5 seconds to 6.5 seconds. 
For each wave period, wave heights ranged from 0.02 foot to 0.18 foot. For 
some of the wave periods the whole range of wave heights could not be tested. 
Enough runs were taken without jet action to establish force as a function of 
wave height curves. Then for each wave height and wave period for 0.6 to 1.0 
seconds three or four different jet-flow rates were tested to establish wave 
attenuation curves. 

All of these tests were run with a water depth of seven inches. This made 


the mooring line scope (ratio of mooring line length to water-depth) equal to 
six. 


Analysis of Data and Results 


The wave and force data were recorded on a six-channel Brush oscillo- 
graph. The incident wave height was calculated by averaging the heights of 
about ten of the incident waves and multiplying the average by the calibration 
factor which correlated distance on the chart paper with actual wave height in 
feet. The transmitted wave height was not so well defined as the incident wave 
height. In most cases it reached a minimum shortly after the jets were start- 
ed. About five of the lowest waves were averaged and multiplied by the cali- 
bration factor to get the transmitted wave height. The wave attenuation was 
then computed using these values and the equation 


Hy - Ht 


= wave attenuation 
; = seaward wave height 
= leeward wave height 


The wave attenuation was then plotted as a function of jet-flow rate. The flow 
rate was determined from the voltage reading with the aid of the calibration 
curve prepared before the tests were started. 

In analyzing the force time records the first two peak forces which oc- 
curred after the wave reached the ship were neglected since they usually were 
due to some starting transient. The length of time for which the run was un- 
affected by side and end conditions was determined and only this length of 
time was considered as the run for the purpose of analysis. The four highest 
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forces during this time interval were averaged and multiplied by the cali- 
bration factor for the force meter to obtain the force in the mooring lines. 
The average of the forces in lines number 1 and 4 (Fig. 3a) were computed 
and plotted as a function of wave height for each different wave period (Fig. 4). 

In order to show the relationship of natural period oscillation to wave am- 
plitude resonance curves for several different wave heights are shown in Fig. 
5. These curves were obtained as a cross-plot from the force as a function 
of wave height curves by taking the slopes of the curves for the different wave 
periods at the various wave heights and plotting them as a function of wave 
period. 

The relationship between wave attenuation and jet-flow rate is given in Fig. 
6 for several ratios of wave length to water depth. A cross-plot from the 
curves given in Fig. 6 is given in Fig. 7 to show the relationship between the 
ratio of wave-length to water-depth and wave attenuation, for several flow 
rates. 


Fig. 8 shows the effect of initial tension on the natural periods of oscillation 
of the ship model. 


Discussion 


The graphs of wave attenuation as a function of jet-flow rate (Fig. 6) show 
that in the range of wave heights tested and within the range of experimental 
error, the height of the waves at a given ratio of wave-length to water-depth 
does not affect the flow rate required to cause a given wave attenuation. The 
general shape of these curves, although not well-defined, due to the inability 
of obtaining data in the range from 0 to 20 gallons per minute flow rates, was 
known from previous work with hydraulic breakwaters.(2,3) The jet-flow rate 
to cause a given attenuation increases as the wave-length to water-depth ratio 
increases (Fig. 7). For the wave periods above 1.0 seconds the jets had no ef- 
fect on the transmitted wave. This shows that the limiting L/d ratio is be- 
tween 6.60 and 8.50 for the range of flow rates available. The attenuation 
curve in Fig. 6c for the 0.60 second period waves (L/d = 3.06) shows more 
scatter than those shown in Figs. 6a and 6b, especially for the condition with 
no jet-flow rate. It has been observed that the waves produced by the wave 
generator at this low period are not good cylindrical waves. This may account 
for some of the scatter of data for this period. 

Curves are shown for the average force in lines number 1 and 4 as a 
function of wave height (Fig. 4). These are the mooring lines on the seaward 
side of the model ship. This average force was selected to use as a repre- 
sentative value because these two forces were generally about equal and high- 
er than those experienced by the other two mooring lines. The forces in the 
mooring lines with and without jets operating were plotted on the same graph 
(Fig. 4). In general the jet action tended to lessen the ship motion and the 
mooring line forces. However, due to the scatter of data it was not possible 
to determine separate curves for the force as a function of wave height with 
and without the jets operating. The scatter of data is believed to be caused by 
the several modes of oscillation which the model experienced for any given 
wave condition, the inability to produce exactly identical waves on any two 
consecutive runs and the two-branch response characteristic exhibited by 
certain non-linear systems. Two wave trains which appeared to be the same 


could be slightly different and cause the model ship to experience two dif- 
ferent modes of motion. 
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Some non-linear spring mass systems have a two branch response curve. (7) 


A description of the response of this type of system may help to describe 
some of the results found in these tests. For low amplitudes of the forcing 
function the mass will oscillate out of phase with the forcing function. As the 
amplitude of the forcing function increases the motion of the mass increases 
until a critical value of the forcing function is reached. When this critical 
amplitude is reached the motion changes phase and “jumps” to a higher am- 
plitude. If the forcing function is then increased still further the motion will 
stay in phase and will continue to increase as the forcing function amplitude 
increases. However, if the forcing function decreases the amplitude of motion 
of the mass will decrease along a different curve and stay in phase with the 
forcing function until another critical point is reached. At this critical point, 
which is at a lower amplitude of the forcing function than the critical point 
mentioned above, the motion of the mass “jumps” discontinuously to a lower 
value and goes out of phase. Then there is one curve for a forcing function 
which starts low and increases and another for a forcing function which starts 
high and decreases. 

The model hydraulic breakwater exhibited characteristics analogous to this 
system for the shorter period waves. In fact at all but the 6.5 second wave 
period condition there was a critical wave height above which the forces in- 
creased very rapidly with increase in wave height. It seemed that up to a 
certain wave height the forces were of the order of the initial tensions, but if 
the wave height was increased by a small amount the force jumped to some 
value very much higher. There was a wave height range in which two con- 
secutive runs with as nearly as possible the same wave conditions produced 
force records which were very different in magnitude. Observation of the 
model motion led to the conclusion that if the motion was such that the moor- 
ing lines did not become completely taut the forces were comparatively low. 
However, once the model ship pulled against the mooring lines hard enough 
to make them taut the motion and the force tended to amplify. Also the lower 
amplitude motion and mooring line force did not always occur for each wave 
period, but tended to be random; whereas the records for higher amplitude 
waves showed peak forces for each wave period. When the forces became 
relatively high, the mode of oscillation experienced by the ship model changed. 
For the low amplitude forces there were two general kinds of ship motion. 
The first can be described as follows: the ship would be displaced ig the lee 
direction and held so that the tension in lines number 1 and 4 was greater than 
the normal initial tension, and the tension in lines number 2 and 3 was less 
than normal. While the ship was held displaced, it would experience a small 
amplitude motion generally a combination of heave and roll which caused 
small amplitude variations in the force in the seaward mooring lines. The 
other type of motion which the model ship experienced for short wave periods 
(less than 1.2 seconds) and low wave heights was when the model ship seemed 
to be excited into a long period oscillation (5 to 7 seconds) by the first few 
waves of the train. Superimposed on the long period motion was a small am- 
plitude rolling or heaving motion. The two types of motion described above 
were generally experienced for the lower wave heights. The critical wave 
height can be seen on the force as a function of wave height curves (Fig. 4). 
That is, when the forces began to increase rapidly with increase in wave 
height, the motion changed. The two types of motion mentioned above were 
relatively small compared to the motion that occurred for the higher waves. 

Above the critical wave height the motion for wave periods less than 1.2 
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seconds was a combination of pitch and yaw. The natural period of pitching 
was 1.34 seconds as determined by the bifilar pendulum method of measuring. 
This method of measurement does not consider the damping due to the water 
or the effect of the mooring lines. The damping due to the water would tend 

to increase this value, but the mooring lines acting as springs would tend to 
decrease the period. Also due to the non-linearities of the system the natural 
period or resonant period would decrease as the wave amplitude increased. 
Since the motion was not pure pitch, but rather a combination of pitch and yaw, 
it may have a resonant period for certain wave heights in the region of the 
wave periods encountered. 

The wave periods between 1.2 and 3.2 seconds and wave heights above the 
critical caused motion which was a combination of sway and yaw. Below the 
critical wave height the motion was primarily sway for these conditions. The 
4.5 and 6.5 second period waves caused sway at all the wave heights tested. 
Even the lowest wave at the 6.5 second period caused the mooring lines to be- 
come taut. 

For some test conditions some of the mooring line forces were noted as 
negative. This is because the force recorded is the force above the initial 
tension in the mooring lines. Since the mooring lines cannot support com- 
pression these negative forces cannot exceed the initial tension. Negative 
forces recorded for runs without jet action were usually on lines 2 and 3 and 
were for the wave periods less than one second. This was caused by the mode 
of motion where the ship was displaced to the lee. Another cause of negative 
forces was the jet action under certain conditions. For large flow rates and 
low wave heights one or two mooring lines would experience forces less than 
the initial tension, because the difference between the flow rates on the two 
sides of the ship model was great enough to cause a net force on the ship and 
hold it displaced. 

A resonance curve is shown in Fig. 5. The ordinate is a measure of the 
rate of change of force with change of wave height. These curves were con- 
structed by taking the slopes of the force vs. wave height curves for a series 
of wave heights. These data then were plotted for constant wave heights as a 
function of wave period. The range of resonant forces increases as the wave 
height increases. 

Fig. 8 shows the dependence of the natural period of the initial tension. 
This plot shows that the natural period of the system decreases with increas- 
ing tension. The natural periods plotted here are for relatively small oscil- 
lations of the ship model. Any motion which pulled the mooring lines taut was 
damped out before one complete cycle could be recorded. Because of the 


damping, it was not possible to check directly the effect of amplitude on the 
natural period. 


CONC LUSIONS 


1. Within the range of wave heights tested and the range of experimental 
error the wave attenuation for a given jet-flow rate and ratio of wave- 
length to water-depth was unaffected by wave height (Fig. 6). 

2. The jet flow rate to give a certain value of wave attenuation increases 
as the ratio of wave-length to water-depth increases. 

3. In general the motions of the ship model were less violent and conse- 


quently the forces were lower when the jets were running than when the 
jets were not running. 
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4. The motions of the ship model were relatively complicated. Small dif- 
ferences in wave conditions sometimes caused large differences in the 

ship model response. 

There was a critical wave height above which the mooring line forces 
increased very rapidly with increases in wave height. This wave height 
decreased with increased wave period. ‘ 
6. The range of resonant wave periods increased as the wave height in- : 
creased. 


7. The natural period of oscillation varied inversely as the initial tension 
in the mooring lines. 


uo 
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APPENDIX I 
List of Symbols 


= wave-length 


d = water-depth 


H; = incident wave height or seaward wave height 


H, = transmitted wave height or leeward wave height 


= wave period 


= mooring line scope 


Q = jet-flow rate 


initial tension in mooring lines 


wave attenuation 


= deep water wave-length 


Ave= average 


APPENDIX II 
Prototype Prediction 


Assume that an area of ocean is to be protected by a Liberty Ship with a 
hydraulic breakwater on board. The allowable wave height in the lee area is 
3 feet. The waves in the area are 10 feet high and the wave period is 9.3 
seconds. The water depth is 50 feet, The length of the Liberty Ship is 454 
feet, and the model used was 5.25 feet. This makes the ratio of length in 
prototype to length in model equal 86.5. From this we get the following scale 
ratios according to Froude Model Laws: 


Length in prototype ae” - 305 


Length in model 


Joet-flow rate prototype 
Jet-flow rate model 


Mooring force prototype 
Moorine, force model 


Time in prototype = \2 
Time in model 


In order to reduce the wave height the necessary amount a wave attenuation 
of 70 per cent is needed. 

The model wave characteristics corresponding to the conditions above are 
T = 1.0, Hj = 0.115 ft. and Hy = 0.035 ft. 

Then the model average mooring line force from Fig. 4c is 2.25 lbs. and 
the jet-flow rate required is 48 gpm or 0.107 CFS from Fig. 6a. Then scaling 
up to the prototype we get: 

Average mooring line force = 2.5 x 648000 = 1,460,000 lbs, 
Jet-flow rate = 0.107 x 69500 = 7380 CFS 
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THE WELLAND CANAL 


W. A. O’Neil! 
(Proc. Paper 1570) 


SYNOPSIS 


The first major step towards the construction of a deep waterway from 
Lake Erie to Montreal, was realized with the completion of the Welland Canal 
in 1932. This canal connects Lake Erie to Lake Ontario, Fig. 1, a difference 
in elevation of 326.5 feet, and provides a channel 25 feet deep and locks to ac- 
commodate vessels up to 735 feet long, so that navigation may bypass 
Niagara Falls. 

The St. Lawrence Seaway Authority was established by the Government of 
Canada for the purpose of constructing, maintaining and operating all Canadi- 
an works in connection with a deep waterway between Montreal and Lake Erie. 
The additional deepening on the Welland Ship Canal to the Seaway Standards 
is therefore being undertaken under the jurisdiction of this Authority. 


History 


The present Welland Canal is the latest in a series of four which have been 
built in the Niagara Peninsula to overcome the obstacle to navigation presented 
by the Niagara River, its Rapids and Niagara Falls. 


First Welland Canal 


The first Welland Canal was built between 1824 and 1829 as a private enter- 
prise by The Welland Canal Company. The builders utilized natural water - 
courses aS much as possible to save expense and labour. The route led up the 
Twelve Mile Creek from Port Dalhousie on Lake Ontario following its waver- 
ing course to the escarpment. Here a deep cut was made across the height of 
land to Chippawa Creek (also known as the Welland River,) at Port Robinson 
where boats descended the creek to the Niagara River and thence to Lake Erie. 
To obtain sufficient water to operate the canal, the summit level was connected by 


Note: Discussion open until August 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1570 is 
part of the copyrighted Journal of the Waterways and Harbors Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. WW 2, March 1958. 

1. Div. Engr., Welland Div., St. Lawrence Seaway Authority, St. Catherines, 
Ontario. 
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a feeder canal to the Grand River at Dunnville. To overcome the difference in 
elevation required 40 wooden locks, each 110 feet long, 22 feet wide with a 
depth of water of 8 feet on the sills. The sailing vessels usually had to be 
towed up the Niagara against the strong wind and current. Ice jams in the 
early spring also held up navigation. An extension of the canal from Port 
Robinson south to Port Colborne eliminating this bottleneck was constructed 
and opened to navigation in 1833, however, the canal was still fed from the 
Grand River. 


Second Welland Canal 


In 1841 the legislature of Upper Canada purchased the canal and decided to 
enlarge it to 9-foot navigation. By increasing the lift, the number of locks 
was reduced to 27, and they were redesigned and built of cut stone, the dimen- 
sions being 150 feet long, 26 1/2 feet wide, and a depth of 9 feet of water on 
the sills. The Port Maitland, Dunnville Branch on the Grand River was built 
at this time. By 1853, it had become evident that the Grand River would not 
suffice as a feeder for a canal with much traffic, and plans were devised for 
a further reconstruction, so that the canal could be fed directly from Lake 
Erie. It was not until 1881 that this was finally accomplished, the water com- 
ing from the lake at Port Colborne. 


Third Welland Canal 


Confederation of the provinces of Upper and Lower Canada in 1867 placed 
the inland water-ways under Federal jurisdiction. A uniform scale of navi- 
gation for the St. Lawrence route and the Welland Canal was recommended, 
with locks 270 feet long, 45 feet wide, and 12 feet of water on the sills, which 
was later changed to 14 feet. The work of enlarging the Welland Canal to 
these dimensions, was started in 1873 and completed 14 years later. 
The course of the third canal utilized much of the previous channel, but 
climbed the escarpment east of the second canal to Allanburg. It then fol- 
lowed the previous route. Its 26 locks were built of cut stone with lifts of 12 | 
to 14 feet. An aqueduct carried the canal over the Chippawa Creek at Welland. 
This improved and enlarged waterway, together with a similar canal sys- 
tem along the St. Lawrence River from Lake Ontario to deep water at Mon- 
treal, met the immediate demands of inland water commerce until the begin- 
ning of the 20th Century, but failed to keep pace with the ever-increasing size | 
of cargo carriers being built on the upper lakes. Transshipments to small Ay 
ships at the Lake Erie entrance to the Welland Canal was costly, and caused = 
traffic congestion on Lake Erie at the junction of 14-foot navigation with 
greater draft accommodation of the Upper Lakes. To remedy this situation, 
the fourth, or Welland Ship Canal was planned. 


Welland Ship Canal 


The route of the Ship Canal was radically changed from the previous canals 
at the Lake Ontario end, so that it runs essentially in a straight line across 
the Peninsula, a distance of nearly 28 miles, Fig. 2. It leaves Lake Ontario at 
a point about 4 miles east of Port Dalhousie, the northern terminus of the 
previous canals, and follows in practically a straight line, due south along the 
valley of the Ten Mile Creek to the foot of the Niagara Escarpment at Thorold. 
This alignment is maintained in the ascent of the escarpment itself, although 
the route of the previous canals is followed from the top of the escarpment to 
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Port Colborne in a general way, all of the sharp bends of the previous canals 
are eliminated. 


Port Weller Harbour 


As no harbour existed on Lake Ontario at the mouth of the Ten Mile Creek, 
now the northern terminus of the Welland Ship Canal, as artificial harbour, 
Port Weller was one of the necessities of construction. By means of a stand- 
ard double track railway, extending from Lake Ontario south for seven miles 
or more along the route of the canal, the surplus excavation of this northern 
division of the work was disposed of in the form of two embankments project- 
ing north from the shoreline 1-1/2 miles into Lake Ontario, to form an arti- 
ficial harbour. At their outer and inner ends these embankments are outlined 
by reinforced concrete cribs, those at the inner end of the harbour providing 
ample dockage space. Rock excavation placed along the exterior of the two 
embankments protects them against storm erosion. 


Locks - General 


To overcome the difference in elevation between the two lakes, 7 concrete 
locks with an average lift of 46.5 feet have been built. The locks are 80 feet 
wide, 859 feet long between gates and with a depth of water of 30 feet over the 
sills. The usable length for navigation is 715 feet and the widest ship which 
has transitted the canal is 75 feet. 

The first seven miles of the Welland Ship Canal south from Lake Ontario 
are banked by the slightly rising lower level of the peninsula, which is de- 
servedly called the Garden of Canada because of its natural rolling beauty and 
intensively cultivated fruit lands. In this section of the canal are the first 
three single locks, Nos. 1, 2, and 3 which elevate navigation 139.5 feet above 
the Lake Ontario level and bring it to the foot of the Niagara Escarpment. 

The climb up the face of the escarpment, by means of equal lifts is ob- 
tained without any deviation from the direct route. Three locks in flight are 
superimposed immediately one above the other so that in a distance of slightly 
over one-half mile, navigation is elevated another 139.5 feet. So that serious 
delays to navigation do not occur, these locks nos. 4, 5 and 6 are twin locks 
which provide separate passage for upbound and downbound navigation, Fig. 3. 

With a short intervening stretch of canal prism, the last of the seven locks, 
which brings navigation to Lake Erie level, is reached in the town of Thorold. 
In this short reach there is sufficient distance for vessels to pass. 

As it is essential to control the summit level of the canal from Thorold to 
Port Colborne at a regulated level of elevation 569.0 there is at Port 
Colborne a supply weir and a guard lock (Lock No. 8) through which vessels 
are locked between the summit level and the fluctuating level of Lake Erie at 
Port Colborne. Lock 8 has a usable length of 1148 feet. A Guard Gate and a 
safety weir are situated just above Lock 7 to retain the summit level in case 
of accident to Lock 7. 


Locks - Filling and Emptying 


The locks in general have been built with concrete entrance walls at both 
ends. The 46.5 foot lift permits the upper gate to be situated on a high con- 
crete breast wall and thereby prevents this gate from being struck by a ship 
proceeding up the lock at the lower pool level. The pond for filling the lock is 
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situated on the east side at its upper end and the waste water passes over a 
weir provided with stoney valves which are opened automatically, should the 
water in the reach rise too high. The waste water flows down through an 
open raceway into the canal prism below the lock. Each lock wall has run- 
ning through it, longitudinally, at floor level, a water culvert 14 feet wide by 
16.5 feet high with a semicircular top. Both culverts lead to an intake situ- 
ated in the east wall, and the inflow of water is controlled by four intake 
valves of the vertical sector type. Each culvert discharges into the canal 
prism below the lower gate, and the outflow from each culvert is controlled 
by two discharge valves also of the vertical sector type. These culverts de- 
liver or take the water to or from a lock chamber through laterals or ports 
3 feet wide by 4 feet high situated in the lock walls at floor level and spaced 
longitudinally at intervals of 30 feet. 

Flight locks No. 4, 5 and 6, the first two of which take the water directly 
from the lock above, are twinned in order that ships may be passed simul- 
taneously in opposite directions through the flight. The flow of the water is 
controlled in a manner similar to that of the other locks, except that the in- 
takes for the culverts in the centre wall are in its upper end, and take the 
water direct from the entrances of the lock chambers, whereas the culverts 
in the east and west walls are fed from an intake in the east wall which takes 
the water from the pondage. 

It requires about 70 acre feet to fill one single lock with a lift of 46.5 feet, 
so it is essential to provide a large pondage area above each lock to avoid 
lowering the canal water level too much when filling a lock. In order to form 
the pondage for Lock No. 6 an earth dam 3300 feet long with a maximum 
height of 75 feet was built. 


Locks - Gates 


The operating gates used at all locks are of the mitering horizontal girder 
type, sheathed on both sides and built entirely of steel, all leaves being 48 


feet long and having girders 5 feet in depth. The lower gates are 82 feet high, 


each leaf weighs almost 490 tons; the upper gates are 35 feet 6 inches high 
each leaf weighing about 190 tons and the gates of Lock 8 and the Guard Gate 
are 44 feet 6 inches high and weigh 235 tons for each leaf. 

Spare gates of each size have been provided, a pair of lower gate leaves 
being stored under water in a special dock below Lock 1. The remainder are 
stored under water in the pondage area east of the shipping channel above 
Lock 1. 

The leaves swing on a hemispherical nickel steel pintle anchored to the 


concrete, and are anchored at the top to two structural steel anchorage frames 


embedded in the concrete. 

The gates are opened and closed by 1-1/4" wire ropes operated by double 
drum winches housed in buildings on the lock wall. Should the gates of Locks 
6, 7 and 8 be carried away, serious damage might be caused to the canal and 
navigation interrupted for considerable time. Double gates have been pro- 
vided at each end of these locks, thus providing extra protection, however the 
real protection is provided by the lock gate fenders. 


Locks - Safety Features 


Certain safeguards have been provided to reduce the chance of damage to 
the canal from misadventure with shipping. All service gates on the locks 
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have rows of mitre safety castings, staggered at their mitre ends, which per- 
mit a gate leaf to be opened about 4 feet before the two leaves lose their 
mitre contact with each other. The essential protection is provided by 3.5- 
inch-diameter wire rope fenders, situated about 70 feet away and carried 
across the lock chambers by booms. The fenders are designed so that should 
they be struck, the boom will be smashed, but the wire rope will be carried 
along by the vessel, the rope paying out against a friction force on each side 
computed to be sufficient to bring a 25,000 ton vessel, travelling at the rate 
of 3.0 m.p.h. to rest in a distance of about 46 feet. Since the canal was put 
into operation, the fenders have been struck 13 times and have provided the 
required protection for the gates. No gates have been damaged by ship move - 
ments in 25 years of operation. 

Long distance electric light signals are situated above and below a lock, 
and regulations provide that a vessel shall not proceed past the “limit of ap- 
proach” sign until the signal changes from red to green. The signals, fenders, 
bridges at locks, gates, valves etc., are all electrically interlocked to give 
proper sequence of operation. This arrangement not only protects the equip- 
ment of the locks, but prevents disaster consequent upon the operation of 
lock gates or valves in other than the predetermined order. Similarly the in- 
troduction of limit switches on all lock and bridge machinery protects the 
equipment from careless operation. 


Locks - Operation 


The gates, valves and fenders of the locks are all operated electrically by 
remote control from control houses placed one at each end of the lock. The 
control boards in the two houses are electrically interlocked with one an- 
other and the movement of the apparatus is indicated, by signal lights. 

Vessels pass through the canal, including the locks, under their own power. 
To control the lock requires a 7 man crew on duty 24 hours a day, 7 days a 
week. The lockmaster is responsible for its complete operation and he super~- 
vises an operator in each of the control houses and 4 linesmen who moor the 
vessel, 

Electrically operated capstans with local controls have been provided to 
handle the ships lines up over the lock walls. 


Lock Maintenance 


To assure a minimum of delay to navigation, an extensive maintenance pro- 
gram is carried out during the 3-1/2-month winter season. In addition, rou- 
tine greasing cleaning and inspections are undertaken, while the canal is 
operating. As with all machinery, proper lubrication is a chief factor in main- 
taining a long life and a low wear rate. There are a number of widely varying 
lubrication programs on the canal. The lock gate operating ropes are re- 
moved each winter and given a thorough lubrication in a pressure vat system, 
consisting of treatment in a cleaning oil, a core lubricant, and a cable dress- 
ing. 

When the canal was built, it was felt that water would provide sufficient 
lubrication for the gate pintles, however, shortly after operation one of the 
pintles seized in its bushing, shearing the dowel pins holding the pintle cast- 
ing. A program was commenced of jacking up the gates, cutting grease 
grooves in the pintles, and installing grease pipes up to the wall coping to per- 
mit regular lubrication during the operating season. Since this system was 
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provided, there has been no trouble until very recently. One of the gates was 
seizing, and on being lifted, it was found that the grease grooves were full of 
bronze filings, which had been gradually scraped off over the years. This 
was due to the grease pipe becoming plugged and, because of inadequate lub- 
rication, the pintle being scored. A program is now in force whereby one or 
two gates each year are jacked up, and new pintles and bushings installed. 

The tainter valves which control the filling and emptying of the locks pro- 
vide a difficult problem to lubrication, owing to their position and construc - 
tion. After the canal had been opened several years, considerable wear was 
noticed at the yoke pins and main bearings of the valves, and a program was 
instituted of installing grease fittings and cutting grease grooves in the appro- 
priate locations in the bushings. This has lowered the rate of wear consider- 
ably, however, each winter a number of valve bushings, and pins are replaced 
as found necessary at the winter inspection. 

Corrosion protection is a very serious matter on the canal, and a continu- 
ous painting program on structures both above and below the water is carried 
out. A bituminous type paint is used on underwater structures, such as gates 
and valves, which are painted during the winter months, and structural black 
and aluminum paints are used on the steel work not subject to immersion in 
water, such as bridges. 

Each year the interiors and exteriors of twelve to 14 lock gates are 
painted, and 35 tainter valves are painted and reconditioned. This program 
continues at a rate sufficient to provide a new coat of paint for each structure 
every five to seven years. 


Locks - Replacement of Worn Elements 


A number of programs are in continual effect, replacing various worn or 
fatigued members. 

Shortly after the canal opened, in 1932, the anchorage I-bars holding the 
tops of the gates were noticed to have fatigue cracks developing at the shear- 
necks, which had been provided in case a ship should strike a lock gate, so 
that the I-bars would fail rather than pull and warp the A-frame embedded in 
the concrete. 

These I-bars were removed and new stronger ones installed without a 
shear-neck. No further trouble has been encountered from this source. 

The anchorage jaw plates which carry the load from the I-bars to the em- 
bedded A-frame, were noticed to have developed slight cracks. These jaw 
plate assemblies are being replaced, the thickness of the plates being in- 
creased from 1/2" to 1-1/4". Where this has been done, no further cracking 
has occurred. 

The stoney valves which regulate the reach levels through the by-pass 
weirs, consist of flat gate valves raised and lowered vertically by means of 
racks and pinions, the thrust against the face being taken by a roller train. 
The roller trains which consisted of steel rollers 4" in diameter and 8" long, 
with journals running in steel side plates, were reconditioned, using several 
methods. One was the installation of cast iron rollers with mild steel jour- 
nals. Another was steel rollers with stainless steel overlay, and a third was 
by installing Meehanite rollers, with journals running in Oilite bearings 
pressed into the steel sideplates. Owing to the difficulty of examining the 
roller trains, an assessment of these different methods will not be available 
for several years. 
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A program was commenced recently of replacing the 3-1/2"' wire rope on 
the wire rope fenders, as it was found to have corroded badly where it passed 
around a fixed snubbing bollard in a semi-circular groove. The cable is in 
good condition where it has been allowed to hang free in the air, but where it 
has been clamped, or otherwise confined, it has rusted. 

Cavitation and corrosion have caused some disintegration of concrete in 
the lock structure. The deteriorated sections are chipped out and the new 
concrete placed as required during the winter maintenance period. 

The white oak timbers which form the mitre sills have been found to de- 
teriorate through time, and replacements are made periodically. 

Electrical maintenance is now largely confined to routine matters following 
corrective measures taken to eliminate severe induction conditions in the 
original installation. 

Interlocking of the various operations was accomplished initially at the 
main circuit voltage of 550 and gave considerable trouble. The interlocking 
circuit was changed to 220 volts which eliminated this condition. 

Cables in ducts, galleries, and undercrossings have been subjected to sur- 
prisingly little corrosion or cold flow of the sheath, but in the case of cables 
in exposed positions where water was able to accumulate in the conduit and 
then freeze, replacement has been necessary. The original tile shelves in- 
stalled through the cable galleries to carry the cables, pulled away from the 
wall due to the weight of the cables, and have been replaced with steel racks 
fastened to the wall. 


Lock - Gate Lifter 


To handle the gates should any be damaged or need taking out for any pur- 
pose, a pontoon gate lifter capable of lifting the large 490-ton gates is kept 
moored below Lock 1. The gate lifter is a self-contained unit, except that it 
has to be towed from place to place by tugs. Trimming of the pontoon is ac- 
complished by shifting water ballast by means of pumps. 


Canal Prism 


The canal is generally 200 feet wide at the bottom and 310 feet at the water 
line, with the banks being sloped on two to one. To prevent damage to the 
banks by wave action, they have been protected by a concrete wash wall 12 
inches thick extending from 5 feet below to 5 feet above water line where ex- 
cavation was done in the dry and by a layer of broken stone protection at all 
other sections. 

The banks have proven unstable at various locations and minor slides have 
occurred. These have been remedied by flattening the slope by dredging, re - 
lieving the top of the bank or placing adequate drains in the banks to reduce 
moisture content. Bank protection and renovations have proven to be one of 3 
the major maintenance items. 


Chippawa Creek 


As in the case of the previous canal, it has been necessary to provide a 
suitable structure to allow Chippawa Creek to cross under the Ship Canal. 
This has been accomplished by constructing a concrete inverted syphon cul- 
vert consisting of six 22 ft. diameter tubes which pass under the canal. 
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Port Colborne Harbour 


The outer harbour is formed by two breakwaters some 7200 feet long 
placed approximately parallel to the shore with a 500 foot opening through 

which navigation passes. A 2000 foot spur extends outwards from the wester- 
ly breakwater providing protection against south westerly gales. 


Bridges 


Twenty bridges, 14 highway and 6 railway span the Welland Ship Canal. 

7 Eleven of these are of the vertical lift type, giving a clearance of 120' from 
the surface of the water. They are counter-balanced by means of heavy con- 
crete weights and chains. Other bridges are of the rolling lift type. Large 
concrete weights counter-balance the weight of the moving span, as the 
bridge rocks upward to an 80° angle. Only two swing bridges are used on 
the canal for the railroads, since they limit the width of the shipping lane. 

A regular lubrication program is carried out on all bridge machinery and 
wire ropes. 

The ropes which hold the counter-weights on the vertical lift bridges have 
been in continuous use since the canal opened, and last winter the ropes at the 
east end of bridge 10, a railroad bridge, were changed as a part of a program 
of replacement. Under the clamps, which gather the cables together, to line 
up for passage over the crown sheave, the wire ropes were corroded. The 
operating ropes which move the bridge up and down are replaced from time 
to time as wear develops. To replace one of these wire ropes, usually takes 
four to five hours. On some of the bridges these ropes lasted 20 years, on 
others they average 8 to 10 years. The difference has been attributed to care 
taken in keeping the tension adjusted properly and in precautions taken to 
keep the wind from slapping them against the tower legs. 

As originally installed, the counter-weight guides were steel castings 
which rubbed on angles back to back, keeping the counter-weight in line, as it 
travelled up and down the tower. Severe wear was noted on both the guide 
castings and the guide rails, and at each of the vertical lift bridges the guide 
castings were machined out, and replaceable bronze gibs installed which has 
reduced the wear considerably. 

The main power supply to the vertical lift bridges is at 550 volts and is 
transmitted to the span via trolley connectors running on three wires up the 
inner vertical face of one tower. The auxiliary control at 110 volts for the 
traffic lights and roadway gates and the telephone line are taken from the 

‘ moving span by two flexible cables attached half way up each tower. The nor- 
mal life of 10 years is often reduced as this cable catches on tower projec- 
tions in a high wind, despite protective guide wires. 


Electrical System 


All valves, gates, fenders, bridges etc., are operated electrically. A 
transmission line extends from the mouth of Port Weller harbour to Port 
Colborne harbour. A Power House belonging to the canal is situated at the 
foot of Lock 4, containing three 5000 H.P. turbine-driven generators, actuated 
by a nominal head of water of 186 feet, conveyed through an 8' -6"' diameter 
penstock from above lock 7 to below Lock 4, and connected to 200,000 volts 
with 550-volt 3-phase power being used for all main motors. Electric power 
at various voltages to suit requirements is used for lighting and heating. 
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Substations for transforming and switching the 22,000-volt power to the 
various stages are situated at locations convenient to the various canal locks, 
bridges etc. 

The main problem to uninterrupted power is caused by tree branches be - 
ing blown into the line during summer storms. In the 25 years of operation 
only a few transformers have failed, and have had to be rewound or emptied, 
dried and reinsulated in spite of a regular inspection and oil testing program. 

In general, the installed capacity of the power house is more than required, 
and hence each generator is run only approximately 1/3 of the time. This al- 
lows ample opportunity for preventative maintenance. 

The turbines are periodically stripped down and the runners built up by 
stainless steel welding where cavitation has eroded them. In addition, the 
generators themselves have been torn down, the cooling ducts cleaned and 
the winding revarnished. The lack of response with the hydraulic governors 
is causing considerable difficulty in operating and changing machines, and 
new hydraulic governors are to be provided, having an electrical sensing cir- 
cuit. This will allow for a more constant speed and steadier 60-cycle fre- 
quency than is attained at present. 


Telephone System 


A complete self-contained automatic telephone system is in use along the 


canal, linking the locks, bridges, substations, power house, and executive 
buildings. 


Navigation 
Canal Traffic 


A comparison of the traffic through the past and present canals is inter- 
esting. During 1831, the second year of operation of the first Welland Canal, 
210,104 bushels of wheat passed through it. It might be noted that there are 
presently several vessels in the Upper Lake Class capable of carrying in one 
trip approximately 600,000 bushels or nearly three times the total amount in 
1831. In 1901, the total tonnage passing through the third canal was 620,200 
tons but by 1914 this had increased to 3,860,000 tons. In 1932 when the 
Welland Ship Canal was opened, the tonnage was 8,535,000 tons and in 1956 it 
had increased almost three fold to 23,066,000, Fig. 4. The ship canal was de- 
signed to allow vessels to make a complete through passage in 8 hours. With 
the present tonnage and number of vessel movements in 1956 the average time 
for an upbound vessel was 8 hours 51 minutes and for a downbound vessel 10 
hours 43 minutes. The larger boats over 640 feet in length require an aver- 
age of 9 hours 14 minutes for an upbound passage and 11 hours 48 minutes 
downbound. The greater length of time for the passage of the larger ships 
has been attributed to the difficulty in safely manoeuvering these vessels in 
the canal, particularly when approaching lift bridges and locks. The larger 
vessels do not respond as quickly as the smaller ones and consequently the 
speed of these ships is reduced by the master so that canalling may be safely 
carried out. 

Cargoes are as varied as the vessels which can be seen in the canal. The 
main items, however, are grain, ore, coal, and petroleum products. Pulp, 
sulphur, limestone and “package freight” are large secondary cargoes. 
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No ship, whatever her home port, pays any toll for the privilege of using 
the canal; but captains must have a “let pass,” obtained either by agreement 
before the navigation season, or at the canal office, to facilitate passage. 

Though tolls were abolished in 1905, a linesmans service fee is charged. 
Ships up to 262 feet in length pay $25.00 while larger vessels pay double that. 
In the old days prior to construction of the Ship Canal, casual labour was em- 
ployed to moor and cast off the lines, the men following along the tow paths. 
Since this system was impossible in the Ship Canal, men were engaged by the 
canal to do this work and the linesman’s fee was established to pay for this 
service. 


Deepening to Seaway Standards 


The present deepening program being undertaken by The St. Lawrence 
Seaway Authority is required to provide navigation facilities on The Welland 
Ship Canal equal to those which will be available upon completion of the 
present construction project, from Lake Ontario to Montreal. This necessi- 
tates the removal of 1,324,000 cu. yd. of common excavation and 1,100,000 cu. 
yd. of rock in the canal reaches, by dredging or dry excavation to provide 27 


feet of water where the bottom is in earth and 30 feet of water where the bot- 
tom is rock. 


Dry Excavation 


Excavation using normal heavy construction equipment is possible for 
three canal reaches between Lock 1 and Lock 4. In these sections where the 
canal can be unwatered during the non-navigation or winter season, three 
separate contracts bounded by the lock structures, were awarded in the fall 
of 1955 for the removal of 735,000 cu. yds. of common excavation. Two of 
these contracts were completed in 1956 and the third is scheduled for com- 
pletion by the spring of 1958. 

In the vicinity of Bridge 4 between Lock 2 and 3 slides in the canal bank 
have developed after deepening. When the canal was originally excavated in 
this area, prior to 1932, severe slides occurred and the top of the bank was 
cut back to relieve the unstable sections. Since that time additional slumping 
of the banks to a minor degree has been observed, however, no interference 
was caused to navigation and the remedial measures consisted of placing 
stone protection to prevent erosion. When deepening was carried out under 
the present program the weight was removed from the toe of the slope and 
shallow slides developed. To prevent further action of this nature, the top of 


the bank is being excavated, thus reducing the surcharge and returning the 
bank to a stable condition. 


Dredging - General 


Deepening of Port Weller harbour north of Lock 1 will be accomplished in 
two stages. The first consists of the removal of silt and loose material by 
suction dredging, disposal being made on the outer face of the harbour break- 
waters. The contract for dredging 340,000 cu. yd. of this material was com- 
pleted in May of this year. 

Dredging of rock shoals which are uncovered after suction dredging will 


comprise the second stage. It is expected that this contract will be awarded 
during the summer of 1957. 
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Another dredging contract was awarded late in 1956 for the deepening of 
3.9 miles of canal southerly from the Guard Gate at Thorold. In this con- 
tract, which is for the removal of 250,000 cu. yd. of common material and 
400,000 cu. yd. of rock, it is extremely important that the turbidity of the 
water in the canal be kept toa minimum. This is a critical matter for the 
paper companies in the immediate vicinity who use unfiltered canal water in 
their process. It could also present a problem for municipalities in the area 
should the turbidity content reach proportions beyond the capacity of their fil- 
ters. To alleviate the condition all common material, including the silt over- 
lying rock, will be removed by a suction dredge. An area must be cleared in 
this manner before drilling and blasting can begin. 

The largest dredging contract which is for 672,000 cu. yd. of shallow rock 
will provide a deepened channel for 4.2 miles of the canal from north of Lock 
8 to Lake Erie. This contract was awarded to two companies in a joint ven- 
ture in the late summer of 1956. 


Underwater Drilling and Blasting 


The rock is being broken by conventional methods using drill boats and 
regular underwater explosives or blasting agents. Each drill boat has four 
spuds and two pneumatic drills mounted on towers which travel on rails along 
the same side of the boat (Fig. 5) thus permitting several holes in a row to be 
drilled before moving the boat. 

An interesting feature is the different methods of blasting used by the con- 
tractors. While similar drill equipment is used, two contractors employ the 
hole -at-a-time method of firing and two use millisecond delays. With the 
first method, individual holes are drilled, loaded with from 8 to 14 lbs. of 
gelatin dynamite and fired. The drill boat remains in place when the charge 
is exploded. When millisecond delays are used, several holes in a row are 
drilled, loaded with a blasting agent “Nitrox,” which consists largely of am- 
monium nitrate, and fired without fear of propagation. Because of the load it ee 
is necessary to move the drill boat at each time of firing. Each of the con- 4 i 
tractors is convinced that his method is equal to or better than that used by 
the other. 

Some difficulty with the delay method has been experienced in that mis- 
fires have occurred. These have been attributed to current leakage in the leg 
wires and extra precautions have been taken to make sure that a ground is not 
made from the wires to the steel scow. In addition the charges are being fired 
using a generator rather than condenser discharge blasting machine so that 
120 volts may be applied to the circuit for the required length of time rather 
than 480 volts for an extremely short duration. This will assure that suffi- 
cient heat will be available to fire the caps. 


Precautions to Prevent Blasting Damage 


Special precautions are required when blasting in the immediate vicinity 
of the lock structures or when near certain sections of the old third canal wall 
which form the present canal wall. Here, close drilling is used to reduce the 
possibility of rock breaking back under the wall thereby undermining it and 
causing a failure. To assure a proper line of breakage one contractor is in- 
stalling in the close drilled holes, sealed metal tubes three inches in diameter 
to provide an intermittent air space in the rock through which the shock cannot 
travel. 
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When blasting rock in water shock pressure waves are set up in the water . 


and can damage nearby structures in contact with it. To obviate such damage 
a method is being employed which was developed by the Hydro-Electric 
Power Commission of Ontario and suggested by Mr. A. Laprairie of Canadian 
Industries Limited who have patented the process. A screen of air bubbles is 
placed between the blast and the structure to be protected by placing a three 
inch diameter perforated pipe horizontally along the bottom of the canal and 
then pumping 15 c.f.m. of air per lineal foot of pipe into the pipe while blast - 
ing. The air bubbles in rising to the surface form the curtain which absorbs s 
a portion of the energy of the blast transmitted through the water. Tests un- 
dertaken by The Hydro-Electric Power Commission in connection with a 
blast set off in the vicinity of one of their power stations at Niagara Falls sug- 
a gest a pressure reduction of at least 70to 1. 
To establish the maximum amount of explosive which may be fired at any 

specific distance from the canal wall, seismic monitoring of the blasts is be- 

ing carried out. Information obtained when blasting at 25 feet from the wall 

in a depth of 30 feet of water indicates that the charge at 50 feet should be 20 

lbs. per delay and at 25 feet 10 lbs. per delay. It is known that existing 

criteria for establishing a safe powder load is not reliable within 50 feet of a 
: structure and therefore trial blasts with small charges at 5 foot intervals will 
be monitored. From the data thus obtained a drill pattern and maximum load 
per delay will be established. 


Underwater Television 


3 Before commencing work in the vicinity of the old third Canal walls at 
| Port Colborne, it was deemed advisable to make an underwater inspection of 
the timber cribs, forming the wall, to ascertain their present condition. The 
survey was planned for a diver, however an underwater television camera 
owned by The National Research Council of Canada was available and it was 
decided to use it to make the survey and at the same time make an appraisal “fe 
of the television system for future work of this nature. The camera in its 
water tight container presented too much resistance to the current to be sat- 
isfactorily handled by lines from the surface and it was necessary to utilize a 
diver to manoeuver the camera into the proper position. A satisfactory pic- 
ture was obtained which was viewed by several engineers and simultaneously 
photographed using a camera focussed on the picture tube. The major disad- 
vantage with the equipment was the extremely small field of view possible 
° with the camera. To obtain a comprehensive picture of a hole undermining ' 
the wall of say 4" x 6' it was necessary to prepare a mosaic of some 12 
photographs. A diver could see the hole and its relation to the remainder of 
the structure from one position. To be used extensively for this type of work : 
. it will be necessary to develop a camera which can provide a clear picture of 7 
G* 


it 


Underwater Excavation 


Standard mucking methods using dipper and clam dredges with dump scows 
are being employed, Fig. 6. In addition one contractor is removing the rock | 
to shore in a unique way. Two flat scows 32' x 170' lashed together side by } 
; ' side have standard land cranes mounted on one end, Fig. 7. The cranes are 
1 equipped with rock grabs which are a special type of clam shell bucket de- 

signed in Germany and operated in the conventional manner. On the land side 
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of these scows, which are placed parallel to the canal bank, is placed a third 
flat scow equipped with a Bailey bridge which acts as a gang plank. The 
bridge can be raised and lowered from the scow. Trucks then travel from 
shore over the bridge and onto the flat scows where they are loaded directly 
by the cranes. The round trip tow to Lake Erie of 38 miles which would be 
required using standard methods is thereby eliminated. 

During the winter, dredging is completed at locations where interference 
would result with navigation during the shipping season. The standard method 
employed is to dredge the rock after blasting with a dipper dredge and scow it 
to a site for rehandling to shore disposal. The scows are dumped and the 
material then picked up from the water by clams and loaded into trucks. Ice 
causes a great deal of trouble while carrying out this operation. In addition 
to tripping the buckets on the dipper dredges it fills the tunnels in one type of 
dump scow making it impossible to open the doors and dump the rock. Asa 
result the scows must be clammed out with a resulting loss in time. The win- 
ter dredging is carried out at a greatly reduced efficiency. 


Navigation During Dredging 


Navigation must continue in the canal while dredging is in progress and 
co-ordination between the contractors operations and the requirements of 
vessels is continually considered. Accurate scheduling of the work along the 
face of docks is of the utmost importance so that the normal operation of 
these facilities will be interfered with as little as possible. As mentioned 
previously, winter dredging is required where interference with navigation 
could be hazardous. This is particularly true at lock entrances where there 
would not be room for a vessel to navigate past the dredging equipment. 


CONCLUSION 


This paper has been presented as a non-technical description and report 


on the main features of the Welland Canal and the work required to deepen it 
to Seaway Standards. 
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ABSTRACT 


Field measurements are presented of forces of up to 40,000 pounds as in- 
duced by two predominant standing waves—those with 1 and 2 1/2 minute 
period which produced water surface slopes of up to 0.0015—in the lines of a 
Ship (USS NORTON SOUND (AVM-1)) of 12,000 long tons as moored alongside 
by unusually taut lines at Port Hueneme, California and near to a node of the 
shorter wave. The two waves acted simultaneously to produce only significant 
surge. It is stated that tight lines tend to reduce the inertial force and that the 
excitation in surge is defined sensibly by the product of the slope of the water 
surface and the weight of the Ship where restoring forces of up to 30,000 
pounds in surge as calculated on this basis are about 1/3 greater than those 
measured. In no case was it considered that resonance was encountered. 


INTRODUCTION 


The reaction of a moored ship subjected in surge to the excitation produced 
by standing waves has been considered systematically at least since Knapp and 
Vanoni (1945) and Carr (1947). The latter likened the moored ship oscillating 
against the restoring force of its lines to that of a mass oscillating with an 
initial velocity against the restoring action of two linear springs with a 
clearance or gap, where the period of oscillation depends not only on the 
spring factor but also on the magnitude of the gap and the initial conditions. 
It is postulated that when the calculated value of this period tends to coincide 
with that of the exciting wave then a condition of resonance is produced with 
relatively high mooring forces. 

Wilson (1951) in a remarkable paper derives an expression for the perio- 
dicity, T,, of the longitudinal oscillation of a moored ship of mass m at reso- 
nance with the exciting standing wave in water of depth d. As rearranged by 


“Note: Discussion open until August 1, 1958. A postponement of this closing date can 
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a. Presented, Oct., 1957, ASCE Convention. 
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Yards and Docks, Dept. of the Navy, Port Hueneme, Calif. 
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the writers, his expression appears to be that for the natural period of oscil- 
lation of a non-linear spring-suspended mass [ Tp = 27 (m/k)1/2 where the 
spring factor k for the lines collectively is equal to D N/4 where D is the 
spring factor for a representative line and N is the total number of lines. 
The k has the unusual dimensions of: force per movement. ) multiplied by 
a factor that expresses the non-linearity of the system and the amplitude, A, 
of the standing wave in the ship’s basin [ (mgA2/kd)¢ where c = (1-n)/2 (1+0).] 
When the lines become significantly taut (n = 1 and k has usual dimensions), 
the second factor becomes equal to unity, the spring factor becomes linear, 
and the general expression of Wilson becomes that for a simple system. He 
showed also the Tp was related to the initial tension in the lines and demon- 
strated that the maximum force in surge to be expected in the moorings as a 
whole was related inversely to this period and directly to the weight of the 
ship and the characteristics of the standing wave, 27 WA/Ty (gd)1/2, which 
expression Deacon, Russel & Palmer (1957) state can be shown to be that for 
the maximum force on the ship if the cables are all tight, the ship is held 
stationary, and the virtual mass of the ship is neglected. 

Wiegel et al (1956) measured the mooring forces induced by short labora- 
tory type “wind” waves on five moored ship models; no particular analytical 
treatment was made of the results. Among other things, they concluded that: 
the importance of initial tension as a factor in the explanation of maximum 
forces had been demonstrated; all of the systems were found to have non- 
linear load-movement characteristics; that in most cases mooring forces 
reach their maximum during the passage of only the second or third wave; 
and that any satisfactory analytical treatment of the problem must include 
coupled motions. 

On the basis that the theory advanced by Wilson (1951) “. . . is so un- 
amemable to analysis that approximations have to be made which spoil its 
many virtues. . .” Deacon et al (1957) present simpler methods for estimating 
the line tensions for the case of a ship with lines that periodically go slack. 
When such a ship is released from rest at an instant when the standing-wave - 
excited water particles are also at rest they indicate by use of simple me- 
chanics that the maximum force required to hold the ship stationary is that 
given by Wilson (1951), and presented herein, multiplied by the factor 1.5. 
Like Wilson, they cite the desirability of tight cables as a means of reducing 
the mooring forces. Joosting (1957) states that while the forces which excite 
a moored ship can be considered sinusoidal, the mooring system generally 
does not have linear extension characteristics, so that no actual resonance can 
be expected, and the general behavior should be examined on the basis of 
forced oscillations. He states further that he has found some anomalies in 
Wilson’s theory when applying it to an important case. On the basis of the 
familiar equation for the forced vibration of a damped spring-suspended mass 
he evolves an expression for the longitudinal movement of a moored ship in 
which the inertia force on the ship is taken to be the product of the mass of the 
ship and the difference between the acceleration of the vater and the accelera- 
tion of the ship. Deacon, et al (1957) have objected to the use of the same rea- 
soning by Wilson (1951). 

An exceptional analytical study of the forces and movements induced by 
ocean waves on a large floating platform moored in about 1200 feet of water in 
the ocean about 300 miles off-shore was made by Continental Copper and Steel 
Industries, Inc. (1956) in which use was made of basic hydrodynamic theory to 
calculate the forces and movement induced by the deep water waves. It 
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appears that a significant restraining force was considered to be present at 
all times so that, as in this study, it was not necessary to cope with the prob- 
lem posed by periodically slack moorings. 


Analytical Considerations 


The condition of force equilibrium on a moored ship in surge, based on the 
damped motion of a simple spring-suspended mass under the action of a 
periodic disturbing force is considered to be: 


inertial - damping - restoring + excitation = 0 (1) 


From a consideration of basic hydrodynamic theory the inertial reaction is 
considered to be made up of two parts: that of the ship and that of the hydro- 
dynamic mass of the ship. This can be expressed as follows: 

9 

inertia force = m (1 + Cy) d*x/ dt2 
where: 

m = static mass of the ship 

Cy = inertial coefficient (ratio of hydrodynamic mass to displaced mass) 


d2x/dt¢ = acceleration of the ship in the surge direction = x 


In most cases, (1 + Cp) is written as Cy so that: 
inertia force = Cy mx 


where Cy is termed a mass factor. 
The damping reaction is based on the well known expression: 


damping force = Cg A P (dx, dt)2 /2 
where: 


Cq = is a drag coefficient, experimentally derived where form is 
considered. 
A = the project underwater area of the Ship for form drag 
wetted area of Ship for frictional drag 
dx/dt = speed of the ship in the surge direction = x 
With cables which are always significantly taut, the restoring reaction con- 


sists significantly of that due to the elasticity of the cables where by the ap- 
plication of Hooke’s Law: (total elongation = P { /AE) 


restoring force = yan tal elong. ( st, A, 


where i represents an item number. 

When there is significant slack in the moorings, as is the usual case, then 
the restoring force is no longer given completely by eq. (4) but rather by one 
of the type given by Wilson (1951): 


* 
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(5) 


restoring force (surge movement) "i 


where: 


k is the spring factor for a single mooring line and n is a numerical 
exponent 
This function must be obtained through a consideration of eq. (4) and the 
characteristics of the moorings in the case of resiliant but significantly 
weightless lines such as ropes. The catenary equation should be used in the 
case of very heavy but significantly non-resiliant moorings such as chains. 

The excitation force from a consideration of basic hydrodynamics is that 
obtained by integrating the proper components of the wave-induced pressure 
over the wetted surface of the ship. It is beyond the scope of this paper to at- 
tempt to.accomplish this analytically; however, as an approximation for waves 
much longer than the ship the David Taylor Model Basin (1954) and Joosting 
(1957) consider the excitation to be that due to a ship sliding up and down a 
sloping water surface such that: 


excitation force = m g dz/dx (6a) 
=m g (27A/L) (sin 27t/T) (6b) 


where dz/dx represents the slope of the water surface in the surge direction 
and A the amplitude of the exciting wave. The above assumes that there is no 
appreciable trim of the ship as it oscillates in surge which assumption seems 
adequate for many cases because the standing waves and combinations thereof 
ordinarily active in harbors are of low amplitude even though very long. 

If the relationships: L = T (g d)1/2 and m g = W are used in equation (6b) 
it becomes: 


excitation force = (27 WA/T (g a)1/2) (sin 27t/T) (6c) 


When t = T/4 or 3T/4—extremes of the oscillations—the sine term becomes 
unity and eq. (6c) reduces to that given by Wilson for the maximum tension in 
the lines as referred to in the Introduction. 

In an environment of standing waves—it is well to consider that moored 
ships will be subject to a maximum variation in water level which is that pro- 
duced by the sum of the various standing waves. Thus, a computation of the 
Slope of the water surface by an expression such as that contained in eq. (6b) 
is not likely to give accurate results, because this equation considers only one 
wave (and that of strictly sinusoidal form) while it is probable that in most 
cases two or more waves will be active. 


Test Facilities and Procedures 


The investigation was conducted in a virtually closed Basin, (Figure 1) 
about 1900 feet long, 280 feet wide and 32 feet deep at MLLW where the tidal 
variation was about four feet. 

The missile launching ship used in the study—the USS NORTON SOUND 
(AVM-1)—is one of about 12,000 long tons (Table I) which has Port Hueneme, 
California as its home port. It was selected for the study because it is 


Figure 1 
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The Harbor and Basin with Ship as moored during the study. 
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Table I Characteristics of the Ship 


Projected area above water line 


26,880 kip pounds 
Mass 


Projected area below water line 


berthed conveniently close to NAVCERELAB; further, it was known that the 
Harbor of the Port contains an interesting system of standing waves which it 
was considered would induce interesting if not high-level mooring forces. 

Because the Basin is part of the Harbor as a whole, where the water in 
both is excited by waves from the Ocean, it is apparent that more than one 
wave system is likely to be active in the Basin at any given time. This is in 
fact the case, because standing waves with periods of about 1, 2 1/2 and 12 1/2 
minutes are known to be active. The Ship was moored alongside Wharf No. 2 
(Figure 1) and near to the node of the wave with one minute period. 

Water-level variations were measured at two points—one near the bow and 
the other near the stern of the Ship (Figure 1)--by means of pressure pickups 
as described by Knapp (1952) placed on the bottom of the Basin. The varia- 
tions due to the long-period waves encountered were reflected directly by the 
output of these pickups, where as it is recognized that this is not the case with 
those of shorter period usually referred to as “wind waves.” The variations 
induced by these long-period waves are important, because as determinats of 
the significant slope of the water surface, they evaluate the major parameter 
in the expression for the Ship excitation (eq. 6a). 

In some cases only the stern water-level pickup was available because of 
malfunctioning of the other; however, on the basis of trends established with 
data from two pickups, it was possible to obtain a reasonable approximation 
to the slope of the water surface from the water-level variation at the single 
location. 

The Ship restraints consisted of seven manila ropes varying in circum- 
ference from 8 to 12 inches (Figure 2) and one spring-lay cable of 1 1/2 
inches diameter. The lines were in various states of condition (Table Il). The 
two surge lines always were kept taut, an unusual circumstance and one which 
served to make easier the task of reducing the force data, because it assured 
that a significant restoring force—later shown to vary almost linearly with 
movement—would be present at all times. This eliminated consideration of 
an unrestrained “run” of the Ship as would have been the case if all the lines 
periodically went slack. 

Mooring forces were measured on a continuous basis by means of dyna- 
mometers of the tube-strain-gage type incorporated directly into the lines 
(Figure 3). Because the incorporation procedure involved a line handling ef- 
fort slightly over and above that required normally, it was generally only on 
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Figure 3. Tension force pickup as used in the Ships lines. 
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weekends, when the Ship did not take to sea, or at other infrequent periods of 
inactivity, that it was possible to activate the dynamometers and collect force 
data. 

Data on the relationship between loads in the mooring lines and their elong- 
ation were obtained by special measurements, some of which were made in 
place (Figure 4). It will be noted that in most cases the relationship is sensi- . 
bly linear where the data presented represents the mean of a number of data 
points which in many cases showed significant scatter due in part to hys- 
teresis effects. 

The outputs of the dynamometers and water-level pickups were transmitted ‘a 
through cables to a photographic-type oscillograph (Miller, Model H) operat- 
ing at a slow paper speed—12 inches per hour—where they were recorded in 
terms of force and water level versus time. 

During one portion of the investigation, measurements were made of the 
on-off forces induced on the Wharf by the Ship, but these proved to be slight 
and further measurements of these forces were abandoned. 

If ship movements were to be measured, it was considered that it would 
have to be on a continuous basis in order to permit selection for analysis of 
only those intervals that proved amenable to analysis. Because a proper pick- 
up of this type for the relatively small movements involved was not available— 
that outlined by Knapp (1952) was found too insensitive, and ordinary survey 
methods were considered too laborious and expensive, direct measurement of 
ship motion was not attempted. Rather, such a determination was made ana- 
lytically on the basis of the calculated elongations of the surge lines. 

No attempt was made to measure mechanically the water currents, because 
no proper meter of the continuous type considered necessary was available. 
Those utilizing propellers and drag spheres proved inadequate. Data from 
floats placed in the area indicated that unidirectional-type currents, including 
those caused by tides, were insignificant. 

Windspeed and direction were obtained by means of a recording selsyn- 
type annemograph which, although located about one-half mile northwest of the 
Ship, was considered to give pertinent data. Wind data was important in the 
negative sense because those intervals during which the wind was at a signifi- 
cant level were eliminated as unsuitable for use in this study. 

The study procedure was to collect data on a continuous basis, then, after a 
review of the records, to select certain significant periods for analysis. There 
were many periods of relatively high disturbance levels when it was not possi- 
ble to collect data because the Ship was at sea. 


Experimental Results 


The experimental results are characterized by an unsteadiness and are 
considered very different from those obtained using uniform, periodic waves 
normal to laboratory-type investigation. Facsimiles of certain oscillograms 
considered typical are presented in Figures 5 through 9. It will be noted that 
the pattern produced by the water-level variations, wherein the major and 
minor peaks shift constantly in a somewhat orderly if not systematic manner, 
seems to indicate the presence of two or more wave trains. The prominent 
long period surge—that of about 12 1/2 minutes—shows considerable variation 
in height but less so in period. Note also that by far the greatest line tensions 
were excited by the shortest of the significant waves—those of about 1 minute 
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Figure 6 Facsimile of oscillogram with waves of 3.50 minute period predominant and a high initial tension. 
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Figure 7 Facsimile of oscillogram with waves of 0.98 minute period predominant and an extra high and extra low initial tension 
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period—because with them the attendent slope of the water surface was great- 
est. Indeed, the two other conspicuous waves—those of about 2 and 12 1/2 
minute period—excited the Ship only a negligible amount in surge, although 
the 2-minute waves induced considerable yaw, as witnessed by the high forces 
present in the stern line (#8). The longer wave induced a heave which at 
times noticeably affected all of the lines. The effect of the initial tension is 
apparent particularly in the Figure 7 data where a quick shift from about 3 to 
22 kips in initial tension in the forward and aft surge lines produced a signifi- 
cant reduction in the restoring force for about the same excitation. 

It is apparent also (Figures 5, 6, 8) that the variations in water level near 
the bow and stern of the Ship due to the standing wave with about 1 minute 
period are 90 degrees out of phase with each other, which is indicative of the 
fact that the center of the Ship is located near to a node of this standing wave. 
(Figure 1-c). This is considered to be a very fortunate circumstance from a 
research viewpoint and influenced to a considerable extent the decision to use 
this Ship for this study. 

The wave spectrum in the Basin was dominated (Figures 5 - 9) by waves 
with periods of approximately 1, 2 1/2 and 12 1/2 minutes. It is not known 
reliably whether or not such long-period waves, assumed to be of the standing 
type, are periodic. It is considered that the various basins in any harbor can 
respond naturally only in a periodic manner depending upon their geometrical 
characteristics; however, it is thought that their excitation is not periodic, be- 
cause it is produced by ocean-type waves. 

In an attempt to bring some order to the results, the period, T,, of the free 
oscillation of the water in pertinent basins of length f was calculated using the 
expression: T = B ¢/n (gd)1/2, where n has a value of 1, 2, 3 etc. for a closed 
basin and 1, 3, 5 etc. for an open basin and B has a value of 2 and 4 for re- 
spectively closed and open basins. These calculations indicate that the funda- 
mental and second harmonic for the Basin oscillation are respectively 1.96 
and 0.98 minutes while the fundamental and third harmonic for the Harbor as 
a whole are respectively 7.50 and 2.50 minutes, the second harmonic being 
supressed entirely. Such calculations do not seem to explain the source of the 
prominent wave with period of about 12 1/2 minutes in Norton Sound slip. As 
speculated by O’Brien (1957), it may be the result of the oscillation of water 
over the very deep submarine canyon which fronts the Harbor. 

No investigation specifically intended to determine the origin of this wave 
is known to have been made. 

In a modest study of the wave environment, a portion of the water-level 
oscillogram obtained from the bottom pressure pickup located near the Stern 
of the Ship was reduced to eliminate the effects of tide and the 12 1/2-minute 
wave. It was then analyzed using the method of residuation as originated by 
Chrystal (1904) and replotted (Figure 10). This procedure involves the off- 
setting of the traces by a certain increment of time to induce cancellation of 
certain components. The original record was residuated first with respect to 
the 2.50 minute wave to expose the 0.98 minute wave and then vice versa, so 
that the result is two separate traces of water level versus time in which the 
waves of 0.98 and 2.50 minute period are respectively predominant. It will be 
noted (Figure 10) that the maximum amplitude of the 2.50 minute wave —0.10 
feet—is about double that of the 0.98 minute wave. From the geometry of the 
Harbor it is apparent that the crest of the 2.50 minute wave as the third har- 
monic of the fundamental oscillation of the water in the Harbor as a whole 
likely is somewhat oblique to the long axis of the Basin; the noticeable 
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tendency of the Ship to yaw when such a wave is predominent may be ex- 
plained on this basis. It is probable that a wave with the fundamental period 

of the Basin—about 1.96 minutes—is also present, but to a minor extent for it 
was not exposed by the residuation technique. This may have been due to the 
tendency toward accumulated errors inherent in the technique. 

The Ship movement was resisted at all times by at least the two surge 
lines—No. 4 and No. 5--which were kept taut such that their tension never be- 
came less than about one kip. This is indicated in Figure 11 which presents 
restoring force—the surge component of the individual line tensions—as a 
function of surge and indicates that movements of about one-half foot result in 
restoring forces of about 20 kips. 

As a standard operation, the tension in the lines for any given tide level 
was varied as desired by the Ship’s crew operating to change the length of the 
aft surge line (No. 5), which operation caused the Ship to move and thus al- 
tered the tension in all the lines. This is illustrated in Figure 11-C where a 
tension of 15 kips in the aft surge line caused the forward bow spring line 
(No. 3) to become active and to operate with the forward surge line (No. 4) 
which was always active. In the case cited there was a decrease in the slack- 
ness of the remaining lines but not enough to induce significant tensions. Only 
in rare cases—as for example during the visit of a personable young lady 
selected by the crew and others as Miss Guided Missile when a sense of gal- 
lantry prompted action to attain maximum ship stability—was the length of any 
line, other than the aft surge, altered. 

When the Ship yawed considerably—as it did frequently when 2 1/2 minute 
waves were predominant—then the stern line (No. 8) was loaded in the order 
of 20 kips. Such intervals were not analyzed because enough difficulty was 
encountered in dealing with significant surge movements without attempting to 
cope with those involving obvious coupling. 

Analysis of the experimental data (Figure 11) shows that the relationship 
between restoring force and Ship movement was continuous and sensibly linear 
in all significant cases, The factor expressing this relationship—designated 
k—was found to vary from 34 to 48 kips per foot as a function of the initial 
tension where a value of about 40 kips per foot seems to cover all cases ade- 
quately. Substitution of this value of k in the familiar expression for the na- 
tural period Tp of the undamped oscillation of a simple spring-suspended 
mass (Tp = 2 7 (m/k)1/2) gives a value for Ty of about 30 seconds. It is obvi- 
ous that this is considerably less than the period of the significant standing 
waves encountered in the Basin during the study, and thus in no instance was 
the condition of resonance encountered. 

An interesting example of the Ship’s behavior is presented in Figure 12, an 
enlargement and elaboration of the basic data presented in facsimile form in 
Figure 5 (point 8). The restoring-force data is the result of summing in the 
surge direction the pertinent components of the measured line tensions. The 
ship movement was deduced from the calculated elongation of the forward 
surge line (No. 5). By graphical differentiation of the ship-movement curve, 
the velocity curve was developed, and from it the acceleration curve was de- 
veloped by the same process. The accelerations were multiplied by the mass 
of the Ship to obtain data for the inertia-force curve; the measured slope of 
the water surface was used in equation (6a) to obtain the data for the 
excitation-force curve. Thus, a sensible approximation of the Ship’s behavior 
in this typical instance is displayed in terms of the equilibrium components as 
given in equaticn (1). 
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Note (Figure 12) that, while the restoring force is always opposite to the 
displacement, the inertial force is in phase as well as the excitation. The 
last circumstance is substantiated not only by the experimental data but also 
by a consideration of the performance of the analogous spring-suspended 
mass when subjected to periodic excitation. It is indicated that, if damping is 
negligible, as can be demonstrated for this study, then there can be no lag be- 
tween the excitation and the displacement. 

In reviewing the effect of damping it will be assumed that viscous damping 
is dominant—that due to wave generation by the Ship is negligible —and that of 
this, the part due to skin friction is negligible since forces of less than one- 
pound are calculated from the expression given by Schlichting (1949) for a flat 
plate in longitudinal flow. This leaves only the drag due to flow separation— 
form drag—to be considered. 

In equation (3) it would appear that the drag coefficient (Cq) should not be 
taken as greater than about 1.0, because for a circular cylinder with axis in 
the flow direction, Rouse (1946) gives a value for the coefficient of 0.90 for a 
condition where the Reynolds Number (N,) is greater than 103 and the length- 
to-diameter ratio is 20. The Ny for this study is about 106 based on the maxi- 
mum particle velocity and the Ship’s beam. The length-to-beam ratio is about 
7. 

The value of x to be used in equation (3) can be estimated on the assumption 
that the Ship is stationary and subject to a standing wave of one-quarter foot 
amplitude (A). In this case the maximum particle speed (Vp) calculated from 
the relationship Vm = A (g/d)1/2 is 0.25 feet per second. From equation (3), 
with the values of Cg and Vp as calculated substituted along with 1150 square 
feet for the underwater projected end area of the Ship, a value is obtained for 
the damping force of less than 100 pounds which is considered negligible. 

At a number of points of maximum acceleration—the extremes of ship 
movement—the mass factor, Cyy, was evaluated using equation (1) and the data 
presented in Figure 12. For accelerations of the order of 0.007 feet per 
second per second, the factor was found to vary between 1.0 and 1.2. The defi- 
nition from the oscillograms of the slope of the water surface—the controlling 
parameter in the expression for the excitation force—was found to be critical 
in the evaluation of Cy. Treatment of a second situation, similar to that 
illustrated in Figure 12, resulted in similar values for Cy. These values 
seemed close enough to unity to permit Cy to be so considered and thus ef- 
fectively to eliminate it as a significant factor. 

The measured maximum values of the restoring forces (as a function of the 
slope of the water surface) are compared in Figure 13 and Table III with maxi- 
mum values of such forces calculated on the basis of equation (1). The mea- 
sured values of the slopes of the water surface were used on the basis of 
readings from both Stern and Bow pickups. In two of the ten cases where the 
Stern pickup only was active, this slope was deduced from the relationship be- 
tween maximum water-level variation near the Stern and the slope as deter- 
mined from data collected when both water level-pickups were active. The 
Ship movement was deduced from a consideration of the elongation of the two 
surge lines on the basis of the particular loads measured in them. 

Values for the maximum acceleration in the Table III calculations were ob- 
tained by assuming the Ship to be oscillating in simple harmonic motion at a 
constant circular frequency, yy, and amplitude, A, such that the maximum ac- 
celeration was equal to AW, That this method of obtaining the accelerationis 
not entirely accurate is indicated by considering the situation which existed at 
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@ = calculated using measured 
water surface slope 
= measured 
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Initial tension = 15 to 22 kips 
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WATER SURFACE SLOPE x 104 


Figure 13 Variation of restoring force with slope of water surface for initial tensions of 
from 15 to 22 kipe. 
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21-04:12:36 Nov. 1954 (point number 8 in Table III and Figure 12) where the 
acceleration as calculated is 0.0057 fpsps, while that obtained by graphical 
differentiation of the ship-movement curve is 0.0075 fpsps. While the differ - 
ence of 30% based on the calculated value of the acceleration is considerable, 
the effect on the restoring force will not be nearly as great—a change of 1.5 
kips in a total of 38.0 kips—because in this case inertia is not predominant. 
It is considered that, notwithstanding the differences cited, the approach pre- 
sented in Table III should be viewed with interest in an engineering field 
where, at present, estimates, especially of the forces induced by waves, to 
within one-hundred per cent or so of the true value frequently are welcomed. 
It will be noted (Figure 13) that the measured values of the restoring forces 
are about three-quarters of those calculated. This circumstance may be due 


in a large part to the liberties taken with the actuai mechanics in the calcula- — 


tion of the acceleration and the equally approximate measures taken to obtain 
the Ship movement. 

When the Ship was tied relatively loosely, as it was in the situations desig - 
nated numbers 1 and 9 in Table II, it was apparent from the analysis of the 
data, if not from intuition, that the ship tended to possess a considerable iner- 
tia force which, when added to that of the excitation, required restoring forces 
that were relatively large—20 to 30 kips. 

It is apparent from a review of Figure 13 data that for any given mooring 
condition the slope of the oscillating water surface is the dominant parameter 
and one with which the restoring force seems to vary almost linearly. 

On the basis of limited data with considerable scatter, an attempt was made 
to demonstrate this further by presenting (Figure 14) restoring force as a 
function of initial tension with the various slopes of the water surface encoun- 
tered drawn as contours in the body of the data. It is apparent that an inverse 
relationship exists for any particular water-surface slope so that it is possi- 
ble to substantiate a conclusion arrived at by Wilson (1951) and Deacon et al 
(1957) among others, and one which may be intuitavely obvious, namely: that 
tight lines tend to reduce the restoring force required. 

As a designer, concerned with the indirect determination of the restoring 
force for non-resonant cases for a ship moored with taut lines in long waves, 
it is speculated that one could cope with the deduction of sensible values for 
the inertial term in equation (1) on the basis of a reasonable estimate of the 
elongation of the ships lines and the amplitude of the water particle movement, 
but that he would be bothered considerably in obtaining, with the same degree 
of confidence, a valve for the excitation. The reason for this lies in the diffi- 


culties inherent in the calculations of the slope of a water surface which likely 
is disturbed by two or more waves. 


CONCLUSIONS 


The following conclusions are drawn from this study of the 12,000 ton Ship 
as berthed with unusually taut lines in a standing wave environment at a loca- 
tion in the Harbor sensibly near to the node of one of two significant standing 
waves, such that the significant movements—those of up to 1 foot induced by 
water surface oscillations at periods of about 1 and 2 1/2 minutes with slopes 
up to 0.0015—were only in surge. 


1. The restoring force for initial line tensions of about 20 kips varied di- 
rectly with the slope of the water surface, being about 10 and 30 kips at slopes 
of respectively .0004 and .0012. 
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Table [J] Calculation of restoring force on basis of measured slope of water surface and ship movements inferred from changes (n tens 


ee line 
Date Point 1 2 3 s 5 t 7 a 9 10 
te year 1954 1954 195% 1954 19% 1954 19s 195s 1955 
dendm 10 Mov 19 Mov 20 Boy 20 Boy 20 Apr Apr 
hour u 18 1 20 2 oT 
ft. not meas 3.90 3.8 3.80 3.65 oot meas not meas 
ft. 3.35 3.590 3.60 3.80 3.45 3.0 & 1.20 
0.40 0.35 0.10 0.40 0.35 0.35 0.9% 
-- #0 00063 0.00098 420.0008 0.00098 0.00085 0.000% 
not meas 3.20 3.65 4.10 3.35 3.50 3.10 wot amas not meas 
ft. 3.55 3.2 3.65 3.95 3.55 3.75 3.55 4.80 1.85 
ft. -- 0.590 0.20 0.20 0.20 0.25 0.45 0.60 - 
height of cecil) 
Bow : tr.-cr. tt not meas 0.@ 0.85 0.3 0.55 0.9 meas oot mas 
Stern: tr.-cr ft 0.20 0.20 0.15 0.35 0.15 0.x 0.25 0.15 
time of oecill. 
trough: B and semis 15:33 15:00 ones 26:00 
crest : Bani & mente 15:57 48:18 06:00 03:30 07:00 14:00 12109 
aiff; T/2 129 124 24 125 
inertial force 
0.017 0.0117 0.0018 0.015 0.0135 0.0L 
at trough 
af 0.0075 0.00b6 9.000% 0.00% 0.0069 0057 0.0019 
force «= 6.3 3.9 o™M 5.5 5.8 5.3 a6 i 
at crest 
Yea .0o97 0.0047 0.0004 0.0077 0.00004 0.0104 0.00% 0.02) 
force « = 61 3.9 os 6.5 5.4 6.7 18.7 2 
@ovement of shi 
6& & 
tension 12.7 15.5 13.0 15.0 15.0 15.0 65 3.5 3.5 2. 
elongation 0.63 oR 0.8 0.14 ,0.2% 
at trough 
tension 23.5 26.3 6.0 25.2 26.5 25.1 v0 $2.0,0.5 26.0 
elongation 1.06 om 1.10 1.10 1.2 1.20,- 
at crest 
tension 2.0 1.5 21.0 5 2.0 Ab 18.5 
elongation 0.06 0.3% 0.97 0.27 0.35 0.06 0.37 
movement ( a ) 
.diff.in elng 
tr. - initial 0.39 0.25 0.40 0.40 1.0 
er. initial 0.57 0.40 0.45 0.64 1.2 
excitation force 
kipe 17.0 26.4 6.5 26.4 22 6 2.6 13.5 
10.8 32.8 13.2 13.2 16.5 64 
restoring force (calc) 
tr.: -Exc. - 27.4 32.7 27.0 6.8 n.9 26.7 26.1 3 26.9 15 
cr.: -Bxc. kipe 25.3 17.1 43.3 19.7 21.9 na 
restoring force (meas) 
— kipe 21.0 17.8 23.8 9.0 20.8 2.6 19.6 are w.0 1.0 
crest 21.0 13.0 25.0 ws ~.7 


Gteined from plet of aeasured slope height of water leve! o#cilietion at the stern 
Nete: 410 ft « distance betwen Botta: Pressure Pickups at Bow ami Stern of Ship 
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2. Restoring force varied inversely with initial tension for any particular 
water surface slope such that with a slope of about .0007 this force was about 
20 and 10 kips for initial tensions of respectively 10 and 25 kips. 


3. Tight lines tended to reduce the restoring force required to resist any 
e particular excitation force since the inertia force was reduced thereby. 


4. Measured values of the restoring force were found to be about three- 
quarters of those calculated on the basis of the slope of the water surface and 
simple harmonic motion. 


5. The excitation can be considered to be represented significantly as the 
product of the slope of the water surface and the weight of the Ship. 


6. The maximum inertial force—the product of mass and acceleration—can 
be calculated to a sensible degree of accuracy by considering the Ship to be 
oscillating in simple harmonic motion. The value of the maximum accelera- 
tion is equal to the product of the amplitude of the Ship motion—deduced from 
the elongation of the lines as calculated on the basis of the forces measured 
in them—and the square of the frequency of the oscillation deduced from that 
of the water surface. 


7. The restoring force can be considered as the sum of the inertia and ex- 
citation forces where damping and the added mass effect are neglected. 


8. A satisfactory estimate for the natural period of oscillation of the 
moored Ship can be obtained by considering this to be analogous to the natural 
period of the undamped oscillation of a spring-suspended mass. A sensible 
value for the spring factor required can be obtained by linearizing the rela- 
tionship between restoring force and Ship movement. 
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SYMBOLS 


amplitude of wave; one-half vertical distance crest to trough 

= amplitude of ship movement 

= projected area of ship 

= cross-sectional area of a mooring line 
B = constant in the expression for natural period of oscillation of 
basins; 2 for closed and 4 for those with one end open 


Ca = coefficient of drag 

Cm = coefficient of mass 

CM = mass factor 
rw c = (1-n)/2(1+n) where n is a numerical exponent used with moor- 
ing ropes 

D = spring factor for representative mooring line 


= depth of water 
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= Young’s modulus of elasticity 

acceleration due to gravity 

item number 

spring factor 

length of a mooring line 

length of a basin 

= length of a wave 

= static mass of a ship 

= total number of mooring lines 

= numerical exponent used with mooring ropes 

= integer (1,2,3 etc) defining harmonic order and nodality of an 

oscillation in a closed basin 
= integer (1,3.5, etc) defining harmonic order of an oscillation 
in an open basin 

= total force in a mooring line 

= time 

= period of oscillation 

= natural period of oscillation of a system 

= maximum speed of water-wave-excited particles 

= weight of a ship 

= horizontal distance along axis parallel to wharf (surge) 

= horizontal distance along axis normal to wharf (sway) 

= vertical distance along axis normal to water surface (heave) 

= angle between mooring line and horizontal 

= angle between horizontal projection of the mooring line and 
the face of the wharf (the x axis) 

= circular frequency = 27 /T 

= 3.1416 etc. 

density of sea water 

speed of ship in surge = x 

acceleration of ship in surge = x 

Reynolds Number 

= (fluid speed) (representative dimension)/kinematic viscosity 
of the fluid 

= sum 
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CONTRACTOR’S PLANNING FOR NAVIGATION LOCKS® 


F. J. Larkin,! M. ASCE 
(Proc. Paper 1572) 


SYNOPSIS 


A contractor’s planning for all jobs starts in the estimating stage. Some 
feature of the work becomes the key to the solution and dictates methods and 
equipment. This paper deals with job planning for navigation locks and is 
principally concerned with the selection of construction plant The con- 
struction of the Markland Locks Project on the Ohio River for the Corps of 
Engineers serves as an interesting illustration. 


INTRODUCTION 


Markland Locks consists of twin lock chambers, one 110' wide by 600' long 
and the other 110' wide by 120Q' long, and is located on the left bank of the 
Ohio River approximately half way between Louisville and Cincinnati. In ad- 
dition to the lock chambers proper, more than 3000 feet of guide and guard 


walls for the approaches are included in this contract. The lift of these locks 
will be 35 feet. 


Planning the Job 


The contractor, on receiving the invitation to bid, is first concerned with 
the location and physical characteristics of the job. He wants to know some- 

thing of the amount of work to be done, which is revealed by an inspection of 

4 the bid quantities. These quantities will have a definite bearing on the se- 
lection of method, and, therefore, the plant. For example, Markland Locks 


Note: Discussion open until August 1, 1958. A postponement of this closing date can be 
obtained by writing to the ASCE Manager of Technical Publications. Paper 1572 is 
part of the copyrighted Journal of the Waterways and Harbors Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. WW 2, March, 1958. 

a. Presented at the American Society of Civil Engineers Annual Meeting, 
October 14-18, 1957, New York, N. Y. 
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involves some 450,000 cubic yards of concrete, 350,000 cubic yards of rock 
excavation, 670,000 cubic yards of common excavation, almost 5 million 
pounds of reinforcing steel, almost 4 million pounds of structural steel and 
castings, a half-million yards of dredging and over 3/4 million yards of fill. 
These quantities, as well as the location of the work tell the planner whether 
or not it is the type of job suited to his organization, and, if he is still inter- 
ested, he starts looking for the key to the method. In planning a job of this 
nature, the engineer makes his plant layout early. (Fig. 1) Excavating 
equipment, of course, has to go to the excavation site, and hauling equipment 
has to be able to transport from that site to the dump. This phase is a study 
in itself, but, at least, its location is set. Not so with the concrete mixing 
plant, which is involved in what is generally the largest single item, dollar- 
wise, on the job. It may be located close to or somewhat remote from the job. 
It involves primarily the transportation of materials both before and after 
mixing process itself. Transportation then may become that key which will 
dictate a major portion of the plant layout. 

In the case of Markland Locks, no rail transportation existed closer than 
12 miles from the job site. A good highway system existed, but the closest 
approved source of crushed stone was over 30 miles away. A contemplated 
pouring rate as high as 1800 cubic yards of concrete per day would have re- 
quired a considerable fleet of trucks to operate and maintain. An alternate 
scheme was the trucking of the aggregates from the nearest rail head. The 
location of the job on the Ohio River, however, suggested moving the materials 
by relatively inexpensive water transportation methods. This scheme is most 
attractive if materials suppliers are conveniently located with respect to the 
river and have barge loading facilities. In this particular case, one supplier 
of sand and small aggregates was located almost on the river bank but had no 
dock. The supplier of large crushed aggregates was located a few miles in- 
land and also had no dock. One cement company who could furnish about half 
of the 400,000 barrels of cement involved was located on the river and had a 
fleet of towboats and barges. The other cement producer was located about 
4 miles inland. In all cases the contractor, being especially equipped for 
river construction, was able to overcome these difficulties, and aggregate 
movement to the job by water was accomplished. Involved were the con- 
struction of an aggregate loading dock near Carrollton, Kentucky, and a simi- 
lar one at Madison, Indiana, as well as a cement loading facility on an existing 
dock at Jeffersonville, Indiana. 

Delivery of materials at the site is only part of the battle. The next steps 
are unloading, storing, processing the aggregate, mixing the concrete and its 
transportation to the forms. With the lock located in the river within a coffer- 
dam, it seemed an unwarranted expense to transport the cement and aggre- 
gates up the hill to the top of the bank and bring them back down again in the 
form of concrete, but such a solution was possible and was considered. The 
primary disadvantage was, of course, the cost of handling materials the extra 
distance. The advantage was relative freedom from high water, which in the 
Ohio Valley can be substantial and of long duration. Of course, if the job is 
flooded out, concrete cannot be placed anyway, so the consideration is simply 
one of keeping dry. Contour maps revealed general top-of-bank elevations 
about Elevation 460, with certain permanent fills going up to Elevation 470 and 
485. A study of hydrographs for the locality revealed one record flood as 
high as Elevation 484, while the duration curve for a 15-year period showed 
the river exceeding Elevation 455, the minimum specified top of the cofferdam, 
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for only 2-1/2% of the time. Storage on top of the bank did not preclude flood- 
ing but protected most of the plant involved. 

A plant located near the material harbor was a logical idea, (Fig. 2) but 
several things had to be resolved before this decision was reached. Usually 
at this point the planner is aware of, or is looking at the plant he has available 
or can buy. In this case, the contractor’s floating plants were of such size as 
to require two of them to suit the capacity required. They were not adapted 
to washing and rescreening the aggregates as required for the crushed rock; 
neither were they capable of passing all the specifications concerning control 
and combined autographic recording of the batched materials, a refinement 
which was specified by the Corps of Engineers. To render them capable of 
passing these specifications required a considerable expenditure of money and 
was particularly difficult to include in the existing plants. 

At this point it is well to consider another element that leads the contractor 
to make certain decisions. On this project four concrete mix alternates were 
offered: 


1. 3" Maximum size stone with all Portland Cement. 

2. 3" Maximum size stone with Portland Cement and one-fourth natural 
cement. 

3. 6" Maximum size stone with all Portland Cement. 


4. 6" Maximum size stone with Portland Cement and one-fourth natural 
cement. 


The difference between 3" and 6" stone represented over 40,000 barrels 
of cement in favor of the 6" stone, or about $140,000 saving in cement costs. 
The difference between natural and Portland cement represented $1 per barrel 
for 100,000 barrels or $100,000. The floating plants were practically limited 
to a 3" size stone and one cement by reason of their hull sizes which are 
governed by lock and bridge clearances. These were, therefore, abandoned 

by the planner in favor of an owned land plant of considerable size which need- 
ed only modernization of its batching apparatus and control to meet the speci- 
fications. 

The apparent solution was to erect the land plant in the cofferdam with the 
shortest aggregate conveyor system possible. One prime advantage realized 
was that short, level concrete haulage was then possible, giving the planner a 
choice between dinkey locomotives and trucks. Storage of aggregate was ar- 
ranged in cells built just inside the cofferdam into which the aggregates could 
be unloaded directly by a whirler crane. (Fig. 3) Draw down of these cells 
in a reclaiming tunnel onto a belt conveyor would bring the aggregates up the 
bank through a rescreening plant and on up to the mixer bins ready for batch- 
ing. Cement could be blown directly from barges to bins located in or along- 
side the batch plant. 

The flood threat was met by placing the batching floor of the plant at Ele- 
vation 471, above any likely high water. This placed the mixing floor at Ele- 
vation 455 or cofferdam elevation. (Fig. 4) The contractor considered the 
possibility and decided to take the risk of having to remove such things as 
mixer motors and control and also the motors and control apparatus in the 
reclaiming tunnel in the event of a flood. From a study of the hydrographs and 
duration curves, this risk was not deemed to be of an excessive nature and 
permitted him to take advantage of the cofferdam location for his plant close 
to the material coming in by water. An 11-foot extension to the bin structure 
was necessary in order to get the plant down to the rock elevation found inside 
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the cofferdam, but with the plant inside the cofferdam and well braced to it, 
this did not prove to be an excessive height of plant. A means had to be pro- 
vided to bring the mixed concrete down to the cofferdam floor, and in this 
fashion the planner arrived at his method. 

The general plan as far as the concrete mixing phase having been decided, 
the remainder of the planning work was considerably simpler. Bank exca- 
vation and the diversion of Stevens Creek (Fig. 5), a stream that ran almost 
at right angles into the cofferdam area, was obviously one of the first items 
of work. Stevens Creek had to be turned downstream to come into the river 
well below the cofferdam. The material, loam, clay and silt, was well suited 
for loading with motorized scrapers and the dump was on the various fills in 
the vicinity. One motorized scraper was available in the contractor’s plant 
and two more had to be purchased. A push dozer in the cut and another on the 
dump, along with a road patrol for the haul road maintenance completed the 
excavation team except for a tractor drawn sheepsfoot roller which would be 
necessary wherever the fills were to be compacted. 

Borings indicated that the lower portions of the cut for Stevens Creek would 
be wet and would not provide good footing for the tractor or rubber-tired 
equipment. For this portion of the work, crawler cranes (Fig. 6) equipped 
with 2-1/2 C.Y. dragline buckets were available to load this material to 22- 
ton end dump trucks for hauling to the fills. At the same time that the bank 
excavation could get underway, work in the river could also start, river stages 
permitting. A considerable portion of the half million cubic yards of dredging 
was located (Fig. 7) in the cofferdam area. Much of the excavation material 
was suitable for coffer fill, so it was essential that the planner figure on the 
dredging proceeding with the cofferdam construction. 

A cellular steel sheet pile cofferdam was judged to be the most adequate 
and economical type cofferdam for this type of job. Here the planner had to 
consider the availability of sheeting. He found nearly enough on one of his 
company’s jobs on Lake Superior. It happened to be in the coffer removal 
stage and so the sheeting could be loaded to barge and towed through the Soo 
down Lake Michigan to the Chicago Canal, down the Illinois River to the 
Mississippi with a left turn at Cairo into the Ohio River and up to the job site. 
(Fig. 8) This situation, of course, was particularly fortunate as it was practi- 
cally impossible to book mill orders for sheet piling at the time, let alone get 
delivery on them. In the bidding for the job it was probable that some con- 
tractors would be forced to consider expedient cofferdams. Quotations, for 
instance, were proffered on a cofferdam constructed from interlocking pre- 
stressed concrete cylinders. This idea was most ingenious but required heavy 
lifting equipment for installation and very considerable cost for removal. 

With these considerations established, the planner figured on using two 
small derrick boats for setting and driving sheetpile. A larger dragline 
derrick boat could dredge the cell sites to scows which in turn would be 
dumped alongside the driven cells, if acceptable for fill, or towed to disposal 
if not. Another large clamshell derrick boat was necessary to fill cells from 
the material dumped from the scows. Enough additional dredging capacity 
would be necessary to remove overburden above the rock inside the coffer- 
dam. This material, or the excess over that needed for coffer fill, was slated 
to be scowed away and would have to come out before the coffer was closed. 
All this work, including the building of the coffer, would have to go at a good 
rate from the end of the Spring rises in the river through the Summer in order 
to be able to pump out the cofferdam in the early Fall. The scheme for placing 
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Fig. 5. Aerial view shows excavation and diversion of Stevens Creek to a 
point downstream of the cofferdam. The stream originally emptied in almost 
the center of the cofferdam area. 
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the mixing plant in the cofferdam meant that erection of it could not start un- 
til the coffer was pumped out. Therefore, the plant would have to be gotten 
ready during the summer and additional items for modernization procured and 
all assembled on the job site ready for a quick start when the cofferdam was 
finally unwatered. Before the coffer could be pumped, it would be necessary 
to cap the cells with concrete to prevent scour in case of over-topping. There 
was enough concrete required for this to justify moving in a mixer boat, and, 
Since this concrete was not specification concrete, river-dug material could 
be used. 

This, then, was the contractor’s plan up through the construction of the 
cofferdam, and the erection of the mixing plant. As soon as the coffer could 
be pumped out, (Fig. 9) clean-up of the bottom would be started, followed 
closely by foundation rock excavation. The planner had noted an area of high 
rock at the lower approach and so planned his cofferdam to extend downstream 
far enough to enclose the bulk of this area. This rock then could come out in 
the dry rather than by submarine methods. 

For the rock excavation the same end dump 22-ton trucks could be used for 
haulage. Loading was to be done by two rock shovels. Line drilling was suit- 
ed to track-mounted drills owned by the contractor, and for shot hole drilling, 
a number of wagon drills were available from the job on Lake Superior that 
furnished the sheetpile. A permanent air plant of 2400 cubic feet per minute 
capacity would be needed, supplemented by 1800 CFM of portable air during 
the drilling period. 

For concrete placing, (Fig. 10) the planner having all his haulage on one 
level, thanks to the plant-in-the-hole idea, selected dinkey locomotives of 
which he owned three at 42-inch gage. Two additional were planned for 
purchase, For placing concrete he planned to use one of his 35-ton gantry 
traveler whirler cranes in the shorter land lock. For the river lock two 50- 
ton similar cranes would be needed which would have to come either from the 
second-hand market or be built new for the job. 

For concreting the upper and lower guard wall extensions which were not 
enclosed in the cofferdam area, the planner figured to set up two small whirl- 
ers, one each on the upper and lower corners of the cofferdam, to rehandle 
concrete in 2 C. Y. buckets hauled from the mixing plant. These buckets (Fig. 
11) would be placed on small barges which could ferry the buckets to the plac- 
ing site. All material handling on these wall extensions would be by the same 
small derrick boats which would assist in the building of the cofferdam. 

This paper has covered the highlights of the preliminary planning that 
started with the estimate for Markland Locks. Once the contractor became 
the successful bidder, planning, especially on details, intensified as dates for 
putting the general plan in effect approached. Of course, the word “planner” 
refers to more than one person. Also it must be stated, as Robert Burns once 
wrote, “The best laid plans of mice and men gang aft a-gley.” A little extra 
high water in the Summer held up coffer work slightly while one mill rolling 
of a barge of heavy section sheet pile was required for the last two cells. 

Also a labor tie-up with the electricians slowed down the plant erection so 
that Spring arrived this year before the plant was ready to run. However, the 
contractor got ahead on rock excavation and form building and by May 15, 
1957 he was ready to roll. (Fig. 12) 
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CONC LUSIONS 


Plant layout and construction methods must be simple, fool-proof and flexi- 
ble. They must also result in economical over-all costs. Simple plant with 
expensive delivered prices on material or complicated and expensive plant 
with cheap materials will not ordinarily result in a successful bid. The trick 
is in selecting the most advantageous method before making a detailed esti- 
mate. Since there usually is not enough time to estimate all possibilities be- 


fore bidding, this phase of a contractor’s planning activities is one of the most 
important to a successful business. 


t 
~ 


March, 1958 


Journal of the 
WATERWAYS AND HARBORS DIVISION 


Proceedings of the American Society of Civil Engineers 


BREAKING WAVE FORCE PREDICTION 


R. L. Wiegel! and R. E. Skjei 
(Proc. Paper 1573) 


ABSTRACT 


Laboratory data of the kinematics of waves breaking on a model beach are 
analyzed and the resulting water particle velocities and accelerations are 
used in the Morison-O’ Brien-Johnson-Schaaf wave force theory to predict the 
force field for piles of three different diameters. Comparisons of results 

with predictions from a modified solitary wave theory are given. 


INTRODUCTION 


In a previous paper (Wiegel and Beebe(5) )* the problem of the design wave 
in shallow water was discussed. A modified solitary wave theory was pre- 
sented and it was shown to predict the velocity profile of the water particles 
under the crest of a breaking wave (in the laboratory) with an average devia- 
tion of zero and a standard deviation of twenty-seven percent provided the 
actual wave crest velocities be known. In addition, it was found that using the 
average ratio of the measured wave crest velocity to the theoretical velocity, 
the velocity profile could be predicted with an average deviation of thirteen 
percent. 
Accelerations were not considered in the previous paper, so a comparison 
was made of the accelerations computed from the “forced fit” solitary wave 
a theory and those obtained by the velocity vector-subtraction method explained 
in a subsequent section of this chapter. These comparisons are shown in 
Table I, and show generally poor correlation. 
It appeared that the modified solitary wave theory should be useful to the 
L design engineer, provided the maximum total moment on a pile occurred 
when a breaker crest passed the pile, because of the good correlation between 


Note: Discussion open until August 1, 1958. A postponement of this closing date can be 

obtained by writing to the ASCE Manager of Technical Publications. Paper 1573 is 
part of the copyrighted Journal of the Waterways and Harbors Division, Proceedings 
of the American Society of Civil Engineers, Vol. 84, No. WW 2, March, 1958. 


*See Bibliography at end of paper. 


1. Associate Research Engr., Inst. of Eng. Research, Univ. of Calif., 
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TABLE I. COMPARISON OF THEORETICAL AND MEASURED 
HORIZONTAL WATER PARTICLE ACCELERATIONS 
IN A BREAKER. 


x/n*= 0.6 cS x/n*= 1.0 


Dist. Meas. Theory | Dist. Meas. Theory | Dist. Meas. Theory 
above accel. accel. | above accel. accel. | above accel. accel. 
bottom bottom bottom 


ft. ftysecs 


5-50 
3-25 
2.19 
1.54 
1.16 
1.03 


3.25 
1.90 
1.28 
1.06 


5.18 
2.57 
1.48 
1.06 


6.18 
3.86 
2.57 
1.87 
1.42 


Wwe 
. 


NR 
. 

OW 
& 


ee 

* 


8.46 
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4.12 
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x = horizontal distance in direction of wave advance as masured 
from the wave crest 


h = water depth below wave trough 


‘ 3 
1573-2 ww 2 
Run 
| | 
22 0.55 0.55 = 7.98 | om 
0.45 -1.1 0.45 2.6 7.98 8.94 
0.35 6.2 0.35 13.9 5.50 0.35 15.8 6.92 
0.25 3.9 0.25 10.8 18 | 0.25 12.7 5.57 
0.15 6.7 0.15 8.3 3.25 | 0.15 19.2 4.57 
0.05 -- 0.05 8.6 2.86 | 0.05 13.6 4.18 
0.25 8.5 0.25 15.7 4.99 0.25 _- 6.34 
0.15 11.0 | O15 15.2 347 0.15 21.8 4.83 
0.05 8.6 0.05 het 2-93 0.05 16.2 425 
0.25 5.3 0.25 19.8 6.50 0.25 7.69 
0.15 11.2 Bot 15.3 3.99 0.15 22.6 Ge 
0.05 10.1 | 0.05 16.6 3.02 0.05 22.0 4.31 
| 
25 | 0.65 263 10.65 — _ 
0.55 0 0.55 9.30 — 
0.45 — 0.45 12.2 6.40 0.45 20.2 7.63 
0.35 0.35 0.35 6,32 
0.25 4.3 0.25 0.25 1.8 5-15 
0.15 6.7 1.13 0.15 | O.15 ll.1 
0.05 -- 1.03 | 0 4.18 ‘ 
26 | 0.85 5457 
0.75 = 4.18 — 
0.65 3.0 3.06 
0.55 2.8 2.19 
0.45 10.3 1.74 ; 
0.35 1.48 
0.15 8.9 1.09 ° 
0.05 1.00 | 
* Note: 
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the modified solitary wave theory and measurements for water particle 
velocities, together with the small effect of water particle accelerations 
under the crest. The only paper which to the authors’ knowledge deals with 
the forces of prototype breaking waves (Morison, Johnson and O’ Brien(1) ) 
and which might indicate the degree to which the above condition would be 
fulfilled makes the statement that the time correlation was not precise 
enough to determine the time comparison between the water surface profile 
and the moment history on the pile. Because of the lack of field data it was 
decided to utilize Iversen’s(2) data to plot water particle velocity and ac- 
celeration fields for breakers, and from these data to compute force fields 
and ultimately the total moments about the bottoms for a series of circular 
cylindrical piles. Since the data presented herein were analyzed, a labora- 
tory study has been made by Hall(1). Hall found that the maximum total 
moment of a breaking wave on a pile occurred nearly at the wave crest, at 


least within the accuracy of one frame of a 16 millimeter movie taken at 32 
frames per second. 


1573-3 


Procedure 


The original data consisted of full-size-scale charts of particle velocities, 
plotted as vectors in the breaker region of a series of laboratory waves 
(Figure 1). From all of Iversen’s data only five runs, at varying wave 
heights, had sufficient data for use in this problem. The vertical distance 
from the bottom to the crest of the breaker, Yp, ranged from about 0.4 foot 
to 1.0 foot. The method by which the original data were compiled has been 
described by Iversen(2). From these charts of velocity vectors, acceleration 
vector charts were derived by the following procedure. Two successive 
charts were superimposed and vector-subtraction triangles drawn (Figure 2). 
From the change of velocity in a particular small area, the acceleration was 
computed by dividing the change of velocity by the time interval between the 
charts. The accelerations thus computed were assigned to the intermediate 
stage between the two charts from which the change of velocity was derived. 
These accelerations were plotted separately on the intermediate wave pro- 
file and analyzed into horizontal and vertical components. The original 
velocity charts were also reproduced, and the velocity vectors also analyzed 
into horizontal and vertical components. Lines of constant horizontal ve- 
locity distribution. Also, constant horizontal acceleration contours were 
plotted on a separate series of selected profiles to determine the horizontal 
acceleration distribution for this particular stage of the breaker. 

The velocities obtained in this manner are probably satisfactory. However, 
the interpretation of the “accelerations” is not quite clear. It is believed that 
they correspond most closely to the total acceleration; that is, they include 
both convective and local accelerations. 

Force charts were obtained for an incipient breaker by the following 
method. An acceleration chart was selected for the desired stage, and the 
velocity charts from which this acceleration chart was derived were again 
superimposed. The horizontal components of velocities over a small area 
were averaged to obtain an average horizontal velocity. This process, re- 
peated over the area of the chart, resulted in an average horizontal velocity 
chart for the same intermediate profile for which accelerations were previ- 
ously computed. This intermediate profile was divided into a one-inch grid 
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and horizontal velocities and horizontal accelerations averaged over the one- 
inch squares. The average values obtained were assigned to the centers of 
the squares. 


The following formulas were used to compute the forces on the pile from 
the above data. 


Fp = 4 Cp p Alulu = Cp pDliu 


and 


Fr = Cyp V du/dt = Cy p md*/4) du/dt 


where Fp = drag force per unit length of pile (Ibs./ft.), Cp = 1.6, o = the 
mass density of the water (2.0 slugs/ft.3), A = projected area of the pile per 
unit length (ft.2/ft.), D = pile diameter (ft.), u = horizontal component of 
water particle velocity (ft./sec.), Fy = is the inertial force per unit length of 
pile (lbs/ft.), Cyy = is the coefficient of mass (2.0), V = the pile volume per 
unit length (ft.3/ft.) and du/dt = the horizontal component of total water 
particle acceleration (ft./sec.2). 

Force-per-unit-length values for the same one-inch squares were aver- 
aged and added algebraically over the grid and the lines of constant force per 
unit length were drawn on separate charts for pile diameters of 0.05 foot, 0.1 
foot and 0.2 foot (Figures 3, 4 and 5). Moments were computed for three 
vertical positions in the breaker: the crest, the maximum-force position, 
and the still water level. These moments were obtained by the equation 


a= i 


where M = total moment about the bottom, F = horizontal component of 
force, 1 = length of vertical section for which F was averaged, and S = height 
of midpoint of the section above bottom (lever arm). Moments were summed 
for the length of the pile to give total moment about the bottom. Moments at 
the crests of solitary waves of the same heights were also calculated, using 
the modified solitary wave formula 


C= 0.81 /g 


and Table I of Reference (5) for the horizontal component of water particle 


velocities, where Y}, is the breaker crest elevation above the bottom. The 
formula 


F= Cop Dlulu 


was used to obtain the force diagram from which the total forces were com- 
puted. The total force has no acceleration component at the crest since the 
solitary wave theory shows zero accelerations along the vertical line through 
the crest. 

For comparison, an alternate procedure was used in addition for one run. 
An original velocity chart was selected, one which had been previously used 
on the second of two charts superimposed for acceleration computations and 
for the average horizontal velocity derivation. The profile was divided into a 
one-inch grid and velocities averaged over the one-inch squares and assigned 
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to the centers of the squares. The advantage here was that original veloci- 
ties were used directly without any averaging inaccuracies being introduced. 
Next, accelerations were computed for the profile by using the change of 
velocities from the selected chart’s preceding chart to the selected chart’s 
succeeding chart. The accelerations were averaged and assigned to the 
centers of the one-inch grid squares as before. Forces were computed, 
averaged and added algebraically, and the total force contours drawn. This 
chart, when compared with the chart drawn by the original method, showed 
very little difference and had an almost identical force contour pattern. 


Results and Discussion 


It is recognized that because the original velocities were derived from 
photographs taken from the side (profile view), the treatment is two-dimen- 
sional and ignores components of velocity which might exist at right angles 
to the plane of the chart. The effect of this particular limitation is to reduce 
the true length of the velocity vectors and to deal only with the chart-plane 
projection of the vector. This two-dimensional treatment will result in 
smaller velocity, or drag-force components of the total force. However, the 
components of acceleration in the planes parallel to the still water surface 
(SWL) are also ignored. The acceleration contour charts show, as a general 
pattern, the zero acceleration line intersecting the wave profile at a point 
slightly above (toward the crest) the intersection of the wave profile and the 
still water level, with negative acceleration (considering direction of wave 
advance as positive) in the area behind on the wave-crest side of this line 
(Figure 6). 

Ignoring this component of acceleration in the negative area results in 
smaller negative calculated acceleration forces and thus result in increased 
total force. As the total force is the algebraic sum of the velocity and ac- 
celeration forces, these two effects tend to balance each other in the critical 
area near the crest of the wave. 

The plotted contours of horizontal components of velocity show an approxi- 
mately vertical zero line near the intersection of the SWL and the wave pro- 
file, with increasing contour lines of positive and negative velocity diverging 
radially from a center near the intersection mentioned. This pattern does not 
change extensively throughout the wave breaking except that an increasing 
positive area develops at the crest of the wave. In these runs a +1.5-ft./sec. 
contour generally extended horizontally from the intersection along the still 
water level, except for Run 26, the largest wave, in which a 2.0-ft./sec. 
velocity contour occupied this position (Figure 7). 

As a wave breaks further, the outer ends of the positive contour curl up- 
wards and a negative area follows. Generally, the maximum positive veloci- 
ties observed were 2 to 4 ft./sec., except in the case of Run 26, the largest 
wave, which had values of from 5 to almost 8 ft./sec., the maximum negative 
velocities observed ranged from 1.0 to <2.0 ft./sec. 

The acceleration contours, as previously mentioned, have a general pat- 
tern in which the zero value contour interesects the wave profile at a point 
slightly above the intersection of the SWL and the wave profile. This zero 
value contour slopes back toward the rear of the wave at varying angles to 
the vertical, ranging from an estimated 40° to about 70°; the lower ends of 
the zero contour generally curved back toward the vertical or past the 
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vertical in some cases (Figure 8). 

Small well-defined areas of maximum positive acceleration occurred 
slightly forward of the intersection of the still water level and the wave pro- 
file in all cases. The maximum positive accelerations observed were values 
of about 20 to 25 ft./sec.2. There did not appear to be any well-defined cor- 
responding small negative acceleration areas. Maximum negative accelera- 
tions observed were generally slightly smaller than the maximum positive, 
but they extended over broader areas. 


Total forces per unit pile length were calculated for these runs using the 
formulas 


Fr = Cm p du/at 


Fo = 4Cpp du 


and using three pile diameters, 0.05 foot, 0.1 foot and 0.2 foot. The values of 
force per unit length, drawn as contours on the wave profile for selected pro- 
file are included as Figures 3, 4 and 5 (in the manner of Schurman(4) ). 
Generally, the same sort of pattern exists for all runs with a maximum oc- 
curring forward of, or at the crest of, the wave. For Run 26, a large maxi- 
mum area occurs at the crest. The original data, in the case of this wave, 
began with the wave well along in the process of breaking. The contours 
drawn for the 0.05-foot pile diameter, and with the same 0.5-lb./ft. interval, 
showed wider spacing than for the 0.1-foot diameter pile, of course, and also 
showed a shift of the center of the maximum force back toward the crest and 
the maximum velocity center, which is to be expected since the drag force 
becomes relatively more important as diameter decreases. The opposite 
effects are seen in the same contour maps, plotted for a 0.2-foot diameter 
pile. The center of maximum force shifts forward toward the maximum- 
acceleration center. Also, the contours are more closely spaced and have a 
considerably more vertical orientation. 

Total moments about the bottom were calculated using the measured hori- 
zontal velocities and accelerations at three locations: at still water level, at 
the approximate location of the maximum total force, and at the crest. These 
were compared with the moments calculated using the procedure described 
in the previous section. The results are presented in Table II. 

The results are quite close except for Run 26, (to be discussed later). It 
is apparent that the total moment which occurs at the crest due to the long 
moment arm is nearly the same as the total moment which occurs at the 
location of the maximum total force for the 0.05-foot, and the 0.1-foot pile. 
The inertia force becomes more important for the 0.2-foot pile since it 
varies directly with the square of the diameter, whereas the drag force 
varies with the first power of the pile diameter. Therefore, the location of 
the zero acceleration contour and the acceleration values are considerably 
more critical and this method of deriving accelerations undoubtedly is not 
sufficiently accurate. 

The calculated values for Run 26 using the modified solitary wave equa- 
tions fell considerably short of the moments derived from the measured 
velocities. It is believed that this is because in Run 26 the wave already 
is past the breaking point, and the measured value of Y}, is not, therefore, 
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TABLE II. TOTAL MOMENTS ABOUT THE BOTTOM FOR THREE SIZES OF PILES 


Total Moments about Bottom 
Run Wave | Still | Wave |Ht. of | We Computed, using |Calc. | Calc. | 
No. |per. |water|;ht. |breaker - | depth|meas. vel. and jusing using 
iT \depth |Hm™ |Hp |acceleratims at |modif. sol. 
\dm wave |pt.of| SWL/sol. wave 


crest |max. | wave | theory 


force) jtheory | (at | 
| |  (aterest) | crest) | 


1.51 2.23 O j 0.07 -Ol 0.08 

1.98 2.23 25 0.04 0.03 

2250 2022 e 0.02 eOl 0.02 

1.00 2.33 0.07 ° 0.10 
0.26 

1.67 2.38 1.14 


1.968 2.23 Ocl: oa 0.08 0.06 
2.22 |0.06 0.03 
2.33 ° 0.13 O. 0.21 


0.25 
10.20 

0.13 
25 0.24 


wad 2.38 . 4.61 


* Measured in the constant depth portion of the wave channel 
*# Using Yp= 1.5 ft. -—- see text 


representative of the true wave. An attempt was made to extrapolate the 
profile of the original wave, from the smooth portion of the profile below the 
breaker and irregular crest of the wave. This extrapolation gave a steep 
profile and a value of Y} of 1.5 feet as compared with the measured Yp of 
1.015 feet. Using the 1.5-foot value, moments were again calculated for 
Run 26 and the resulting values agreed quite well with those computed from 
the measured velocities, as seen in the table of moments. It is important to 
note that the field measurements in Reference (2) showed this same thing; 
that is, the forces exerted on the pile by a “foam line” of a certain height 
were considerably higher than the forces for a breaking wave of the same 
height. 

Recent laboratory measurements by Hall(1) show that moments given in 
Table II are fairly reliable, but the scatter of the measurements indicate that 
it would be better to use the solitary wave velocity rather than 0.81 times the 


solitary wave velocity in computing the water particle velocities to be used 
in the wave force equation. 
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CONCLUSIONS 


It appears that the solitary wave theory should be useful for the prediction 
of total moments about the bottom for piles of small size in breaking waves. 
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CELLULAR COFFERDAMS AND DOCKS®* 


Discussion by Samuel Heyman and j 
Cevdet Z. Erzen 


SAMUEL HEYMAN,* J.M. ASCE.—Mr. Cummings is to be complimented 
for presenting some new ideas on a difficult subject and for publishing 
the results of his tests. Additional tests on cellular cofferdams are re- ; 
quired for determining where in the fill the actual plane of failure occurs; : 
however, Cummings’ results are a valuable addition to the available data on 
the behavior of cellular cofferdams. 

In his paper, Cummings states that ¢ was determined from the force re- 
quired to make the model slide. After making several tests, he divided the 
average horizontal force required to make the cell start to slide by the weight 
of the cell; this yielded the tangent of the angle of internal friction, ¢, of the 
soil. After several trials, he decided that the average angle of internal fric- 
tion was 38° 56'. Cummings uses this ¢ value of 38° 56' to compute the 
theoretical strength of the test cell and then shows how well the actual test 
results agree with the theoretical. This is unjustified since the angle of in- 
ternal friction was determined from the test results themselves and not by an 
independent method. 

The following analysis has been made to demonstrate two ideas. The first 
is that all the major methods of cellular cofferdam design become invalid at 
extreme, unrealistic values of ¢. The second is a comparison of the various 
design methods in the range of practical values. 

To demonstrate the first idea, the resisting moment of the cohesionless 
fill as determined by Cummings has been compared with the resisting moment 
of the cohesionless fill as determined by Terzaghi. (1) Terzaghi states that 
the unit shearing resistance along a vertical plane of the cohesionless fill 


is equal to ; Yh2 C tang where C is the horizontal earth pressure coeffi- 


cient. Terzaghi analyzed this force, and from the dimensions of cells which 
have worked satisfactorily decided that C “must be considerably greater than 
the active Rankine pressure” coefficient and “for the middle part of the fill in 
the cells ranges between 0.4 and 0.5.” The author, however, in the discussion 
of Terzaghi’s method considers it as using the active Rankine pressus along 
the vertical plane. In the published discussions of Terzaghi’s paper, (2) 
Professor Krynine demonstrated (and Terzaghi sencenedl®) by a Mohr’s 
Circle analysis that this coefficient C should be equal to 


a. Proc. Paper 1366, September, 1957, by E. M. Cummings. 
1. Head, Sails Dept., Frederic R. Harris, Inc. New York, N. Y. 
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cos” g 
2 - cos* g 


Terzaghi’s formula for the moment of resistance of the cohesionless earth 
fill in the cofferdam against vertical shear is 


Ms= C tan 
When the relationship of b = 0.85h is assumed, 


Ms = 0.283 ¥ C tan 


For any given cell 0.2837h3 isa constant; and C tan ¢, and therefore Mg, 
varies solely with ¢. 


The moment of resistance of the sand as expressed by Cummings in 
formula 6 is equal to 


Substituting c = b tan ¢, a = h-b tan ¢, and Terzaghi’s relationship 
b = 0.85h, this formula reduces to: 


Ms = (0.283 ¥ h?) [ 0.362 tan (3.53-tan tan = 0.283% h° M tan 


Where: C=M= 0.362 tan (3.53-tan #) 


“M” will be called the Cummings Coefficient. Since both Terzaghi’s and 
Cummings’ expression of Ms can be reduced to the same form, it is suffi- 
cient to compare the variation of the “C tan ¢” term in each case. 

On the following graphs is a comparison of the various “C” coefficients 


that are used in the design of cofferdams. The following designations are 
used on the graph: 


C=M= 0.362 tan (3.53-tan f), 50-4, 


2-cos 
and 


C=astan®(45- 


The term C tan ¢ rather than C alone has been plotted against ¢ because C 
tan @ is the entire parameter which varies with ¢, and because it is consist- 
ent with the author’s C tan ¢ presentation. As can be determined from this 
graph when ¢ = 90°, the Cummings’ Coefficient “M” times tan ¢ approaches 
minus infinity. Describing Terzaghi’s method the author, believing that 
Terzaghi’s method utilizes active Rankine pressure for failure along a verti- 
cal plane through the fill, points out that when ¢ becomes 90° “the shear, 
and thus the resisting moment against tilting, would be zero, as though the 
cells were left filled with water.” Using the Cummings’ method, when 


4 

2 3 
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¢ = 90° the shearing resistance of the material within the cell would be 
infinitely worse than plain water. 
As stated above, the original formulae were not set up for such extreme 
values of ¢. In the case of Cummings’ formulae the values of “a” and “c” 
can never be greater than “h,” which would limit the value of @ to approxi- 
e mately 47° using the b = 0.85h relationship. This demonstrates that all 
theories give impossible results for impossible @ values. 
The values of C tan ¢@ in the range of practical values of @ for a cohesion- 
less soil are shown on the accompanying charts. As can be seen, the 
. Cummings’ design method is the least conservative for @ greater than 30°, 
and that using the coefficient of active pressure times the tangent of ¢ for 
shear along a vertical plane is the most conservative design method by a 
wide margin. For @ less than 30° Krynine’s method is less conservative than 
the Cummings’, and at ¢@ = 30° the maximum of Terzaghi’s range is within 
20% of both Cummings’ and Krynine’s values. Because the differences in 
design coefficients, especially near @= 30°, are smaller than the normal 
factor of safety used in design, cellular cofferdams have been successfully 
designed by the various methods discussed herein. 


REFERENCES 


. “Stability and Stiffness of Cellular Cofferdams” by Karl Terzaghi, 
Transactions ASCE Vol. 110 1945 P. 1083. 


. Ibid. Discussion by D. P. Krynine P. 1175. 
. Ibid. Discussion by Karl Terzaghi P. 1187. 


CEVDET Z. ERZEN*—Mr. Cummings considers horizontal planes in a 
cellular cofferdam as the planes of failure due to shear in soil and introduces 
a theory based on this assumption. This concept seems justifiable when a 
model with 3/4 inch thick vertical planks which correspond to very stiff piles 
in the prototype is observed. Inasmuch as the back of the model stays 
straight, thus producing deflection varying linearly, considering the theories 
of Rankine and Coulomb there is no reason why the planes of failure in the 
cell must be horizontal. 

The use of 3/4 inch planks in 19-1/8"' x 24-1/2"' x 24"' cells which cor- 
responds to very stiff piles in the prototype creates doubts as to the exact 
duplication of the behavior of cells with flexible sheet piles. 

Dimensional analysis of such a system renders 


F f,_F B 


EL EL?) 


Us 


u = deflection of sheet piling 
F = applied external load 


B= flexural rigidity of pile i 
h L = length of pile | 
E = modulus of elasticity of material | 
*Frederic R. Harris, Inc., New York, N. Y. | 


| 
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Denoting the scale factor of the model by K, there are obtained 
KF 


Taking the modulus of elasticity of white pine at 1,250,000 psi, the follow- 
ing are found: 


30,000 _ 
Kp = = 24 
Kp = 24 Ky 2 
Kp = 24 K,4 
Ky = Ky, 


Assuming that the cofferdam is ten times the size of the model 
Ky = 2,400 
Kp = 240,000 
K, = 10 

In other words, the line load applied on the cofferdam is 


F = ee (2,400) = 240 kpf 


The flexural rigidity of the sheet pile is found to be 
3 
= (240,000) i2 (12) (.75)* Em 


where the subscripts p and m denote prototype and model respectively. 
Solving for Ip there is found: 


Ip = s0/100 = 4,220 in 4 per ft. 


This is equivalent to using 24 WF160 instead of the MP101 or MP102 sec- 
tions which possess almost negligible moment of inertia. 

Since the bending of such a heavy section is insignificant, deflection varies 
linearly from zero at the bottom to a maximum at top. If it were assumed 
that this would produce slip planes due to shear on horizontal planes, then the 
theory introduced in the paper would be correct. However, according to 
Rankine or Coulomb theory the planes of failure are different from the 
horizontal planes. 

It is further attempted to show the existence of two zones, namely, active 
and passive, by means of pull-out tests on rods embedded in sand. In the 
opinion of this writer, the proximity of the outer rods to the boundary must 
be considered in evaluating the results, especially in tests indicated in 
Figs. 6 (b) and 7 (b) where no external force is applied. Neither the plastic 
nor the elastic concept of stress in circular discs can justify the substantial 
difference in loads in Fig. 6 (b). Therefore, it seems logical to consider 
the effect of the boundary on the stresses acting on the rods. 


F 
ad 
Vat 
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In the paper, the resisting moment of the cell against overturning forces 
is considered to be the sum of the resistance due to interlock friction and 
that due to shearing strength of the soil. Clearly, the moment of resistance 
of each pile is PLfb, if there exists no stress on the cross sectional area of 
diaphram or if this stress is constant. However, there exists bending 
ae stresses on the cross sections of the piles even when the cells are not em- 
bedded (as in the case of cells on rock), by virtue of the bending stresses in 
the soil which are transferred to the diaphragms through shear acting be- 
tween the piles and the soil inside the cell. Therefore, the moment of 
. resistance is not due to interlock friction acting at the ends of the diaphram, 
but it is due to the vertical shear at center which is some portion of the 
triangular area under the cell arising from the bending of the cell. Thus, 
this moment is 


2 
Mf = 3 Pfb (1) 
The difference between the external moment and the above value must be 
resisted by shear in the soil acting on a vertical plane. 

The writer has used the concept of the existence of plastic and elastic 
zones within a cell in the design of cofferdams. Such an analysis should be 
consistent with the laws of statics and Coulomb’s theory. Irrespective of the 
dimensions of a cell there may exist an elastic zone. The forces in the plas- 
tic zone are given according to Rankine’s or Coulomb’s theory. Once these 
forces and the two zones are determined, the stresses in the elastic zone 
may be calculated by treating this section separately with all the known loads 
from the plastic zone acting on it. 

The failure of the cell is possible if the maximum shearing strength of the 
soil in the cell is exceeded unless the overturning moment is less than the 
moment given by eq. (1) and also the sliding along its base is prevented. 
These principles may be explained by considering a cell resting on rock and 
holding back a fill. 

In Fig. 1 (a) is shown such a cell. The area ABCDA is in a state of plastic 
equilibrium where the shear in the vertical and horizontal planes is zero. 
The active pressure on plane AB is 


Pa = y gh?Ka 

where subscript 2 refers to the plastic zone and similarly subscript 1 will 

¢ refer to the fill outside the cell. This pressure is transmitted through the 
diaphragms to the back piles of the cell and is resisted by the elastic zone. 
If the elastic zone is ADFA, the total force acting on DF is 


1 2 2 

3 Ka + Y 9 (1-c,) 
and the shear on the base AF is 
1 
2 
If G is the mid point between A and F, and assuming the shear is equally 
distributed between AG and GF, the total normal force on plane GJ is found 


from a free body diagram of forces acting on ICFGI or IGABI. Summing up 
forces in the X direction there is found: 


S=5 7, h?Ka (2) 


7 
4 {| 
7 
7 
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2x, 2 2, _S 
+5 Yoh Ka (1 -a,) -3 (3) 


where the second term is the force due to the pull in the diaphragm between 
G and J. The maximum shearing stress in the soil may be found from Mohr’s 
circle of stress. For this purpose, normal and shearing stresses at mid _ 
point between G and J must be shown. The average normal stress on GJ is 1 
N 


and the normal stress on a horizontal plane is 


h + ah 


The shear on plane GJ is found from the moment of the fill about the base. 
Since the external moment is resisted by the soil and the sheet piling, the 
moment due to soil is given by 


Ms = M - My 
and the vertical shear on plane GJ is 
Ms 


The average shear on this plane is 


GJ 


Knowing 0x, oy and 7, one may obtain the maximum shear stress and 
the normal stress on the plane where this shear exists. These stresses may 
be also found from the formulae 


T max = (A 2, 


on = 


If T max>on tan @¢, then failure occurs internally and Ms becomes zero. In 
this case, in order that the external forces of the fill may be still resisted, 
the moment of the interlock friction must be greater than the external » 
moment. 

The boundaries of the elastic zone are determined from a consideration 
of the laws of statics. The elastic zone is not necessarily a triangle ADFA 
as shown in Fig. 1, but it may have the boundary AJEFA. To find this zone 
a free body diagram is drawn of a vertical section through H as in Fig. 1 (b) 
where the volume of soil HEFKH is considered separated from the sheet piles. 
Therefore, only the forces acting on the soil above HE and those acting on 
the sheet piles are considered. 


Summing up the forces in the X direction there is found as the first 
condition 


2 2 2 2 
Pay + Pag (Of - OH)>Pp (aE? - 
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and the second condition is that due to overturning. Taking moment about E,: 


h 
+ Pag (ak? - aH? ) y> 


P, (OE + Pagf (1-p°) (5) 


where the last term is the friction force in the interlock. These are the 
two equations which determine the elastic zone. 


. As an example consider the cell shown in Fig. 7. With Y = 91.9 pcf, h=2', 
¢ = 38° - 56', the following are found: Tan ¢ = 0.81, Ka = 0.228, Kp - 4.39, 
Py = 803 lb., Pap = 42 Ib. 
Taking ap = 1 there is found oy” = 0.428. The first terms in eq. (5) is 
the moment of the external force about E, or according to Fig. 7 


Ry > 2 (P, - Pag) (aE2 - aH?) 


= 2x 761 x 0.572 x .321 = 279'# 


which is the maximum moment that can be developed. The maximum value 
of R is 485# which is determined from the consideration of sliding on the 
base. Maximum moment of 279 ft. lb. is a good approximation to the actual 
test results and the curve Ry = 279 ft. lb. follows the curve obtained from 
tests very closely. 

Shear on the bottom is R = 200 lb. Assuming this to vary according toa 
parabola, the shearing force within a distance of 3 inches (where the rod is 
located) is 13 lb. Force on the back of 2 ft. wide cell is R + 2Pag or 284 lb. 
The normal force at three positions of the rod are 


Ny = 284 - 13 = 271 lb. 

No = 284 - 100 = 184 lb. 

Ng = R + Pag (1 - 0.752) - 187 = 50 Ib. 
from which the force necessary to pull each rod becomes 


271 
Pi 2-5x 2x 0.81 = 45 lb. 


Po = 30 Ib., Pg = 10 Ib. 


f< Example given in Fig. 12 is treated here for demonstration. The elastic 
zone starts at A and makes an angle of 27.5° with the vertical and changes to 
30° as shown in Fig. 2. First, eq. (4) is used to find if there is any possibili- 
ty of a failure plane through L. Here, L corresponds to point E, and B to H 

: of Fig. 1. In eq. (4), 


Pay ap? = 96 
since Pay aE? + Pay (aE? - aH?) P, (aE? - alt”) 


failure does not occur along BL. Similarly eqs. (4) and (5) are not satisfied 
for failure along ME. Therefore, the elastic zone is bounded by line ABC. 


q 
7 
1 
- 
oH?) = 426 
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Total shear along NM is S = 147 - 75 = 72K and the normal force on KG, 
with 


72 


N = 147 + 87 x .85 - —5 = 185K a 


Therefore 


185 
ox ae = 2.5 ksf 


The vertical pressure mid way between K and G is 
oy = 5.2 ksf 


and 


Ms = 6098 - 623-2 1306 = 4605 


128 
= 4 = 1.74 ksf 


Using these values the max shear is found to be 


Tmax = 2.21 ksf 
and 


on = 3.8 ksf 


Tall. = 3.85 x 0.70 = 2.70 ksf 


since tT all > 7 max. the soil is in a state of elastic equilibrium. However, 
in other cases the maximum shearing stress in the lower zone may exceed 
the maximum allowable shear in which case there exists a state of plastic 
equilibrium through the cell. 

The theory given by the writer determines these zones through the use of 
eqs. (4) and (5) and also establishes the state of equilibrium of the lower 
zone - whether elastic or plastic. If the lower zone is in a plastic state, as 
a first condition, the external moment must be resisted by the moment of the 
interlock friction. The second condition of equilibrium is that due to sliding. 
For equilibrium the sliding forces must be less than the sum of the resisting 
forces and the shear on the bottom due to overburden. Finally, the equilib- 
rium of the overturning moment of the external forces must be satisfied. 


| | 
_3 Ms _ k ; 
v=5 4° = 128 
giving 
| 
a 
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¢. Discussion of several papers, grouped by divisions. 
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